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OJJi Abstract. The generating series of a radix-rational sequence is a rational 

3 ' formal power series from formal language theory viewed through a fixed radix 
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precision of the asymptotic expansion depends on the joint spectral radius of 
Cn ' the linear representation of the sequence; the coefficients are obtained through 

some dilation equations. The proofs are based on elementary linear algebra. 
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1. Introduction 

Radix-rational sequences come to light in various domains of knowledge and 
almost every example has been first studied for itself and by elementary methods. 
As a resu lt there does not exi s t a general theorem about t heir asymptotic behaviour. 
Flajolet ( Flajolet and Golinl . Il994t iFlajolet et al.l . 119941 ) has developed a method 
based on the analytic theory of numbers. The used frame of divide-and-conquer 
recurrences is wider than ours and permits to deal with many examples, but even 
if the underlying idea has a geometric evidence the application remains delicate. 

Yet these sequences are a mere generalization of classical rational sequences, 
that is sequences which satisfy a linear homogeneous recurrence with constant co- 
efficients. In the same manner they satisfy a linear homogeneous recurrence with 
constant coefficients but where the shift n i— > n -I- 1 is replaced by the pair of scaling 
transformations n i— > 2n and n i~* 2n+l (or more generally n i— > Bn + r, < r < B , 
if the used radix is B > 2). The first example which comes in mind ( TroUopd . 



19681 : iDelangd . Il975f ) is the binary sum-of-digits function, that is the number S2(?^) 
of I's in the binary expansion of an integer n. (I Prmota and Gajdosikl ( 19981 ) give 
an extended bibfiography.) The sequence satisfies S2(0) = 0, S2(2n) = S2{n), 
S2{2n + 1) = S2{n) + 1. It may seem to the reader that such a sequence is of limited 
interest, but it appears in many problems like the study of th e maxim um of a deter- 
minant of a n X n matrix with entries ±1 (Clem ents and Lindstromlll965f) or in a 
merging process occurring in graph theory (jMcIlroyl . 119741 ). This example have been 
greatly ge neralized with the number of occurren ces of some pattern i n the binary 
code ( Bovd et al.l . Il989r ). or in the Gray code: Tlaiolet and Ramshawi ( 19801 ) study 
the average case of Batcher's odd-even merge by using the sum-of-digits function of 
the Gray code. Among the sequences directly related to a numeration system the 
Thue-Morse sequence which writes u{n) = ( — 1)^2(") is certainly t he on e which has 
caused the greatest number of public ations (Allouche and Shallit ^ l999[ ). There ex- 
ist varian ts with some subsequences (JNewmanl . Il969t ICoqueti . Il983l ) or with binary 
patterns OBovd et al.l . |l989) other than the simple pattern 1: the Rudin-Shapiro 
sequence associated to the pattern 11 was initially designed to mini mize the L°° 
norm of a sequ ence of trigonometric polynomials with coefficients ±1 ( Rudinl . ll959t 
Shapirol . Il951[ ). The study of the complexity of algorithms is another source of 
radix-rational sequences. The cost c„ of computing the n-th power of a matrix 
by binary powering satisfies cq = 0, C2n = Cn + 1, C2n+^ = Cn + 2. The idea of 
binary powering has been re-employed bv lMorain and Olivoa ( 19901 ) in the context 
of compu tations on an elliptic curve w here the subtraction has the same cost than 
addition. Supowit and Reingoldl (|l983 h have used the divide-and-conquer strategy 
to provide heuristics for the problem of Euclidean matching and this leads them to 
a 4-rational sequence in the worst case. The theory of numbers is another domain 
which provides examp les like the number of odd binomial coefficients in row n of 
Pascal's triangle (iStolarskvl . Il977l ) or the number of integers w hich are a sum of 
three squares in th e first n integers (Osbaldestin arid Shiul . 119891) . 

It is the merit of lAUouche and Shallid ( 19921 . 120031 ) to have put all these scattered 
examples in a common framework. A sequence (u„) is rational with respect to the 
radix B if all the sequences obtained from it by applying the scaling transforms 
n <—>■ Bn + r, < r < B, remain in a finite dimensional vector space V. This leads 
to the idea of a linear representation with matrices of size d if y is d-dimensional. 
For the examples above the dimension d is usually small, say from 1 to 4, but there 
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exist examples where d is larger, like in the work of ( Cassaignel 119931 ) which uses 



rf = 30. In such cases a general method of study is necessary. 

Our aim is to provide asymptotic expansions for radix-rational sequences. Be- 
cause the sequences under consideration may have a chaotic behaviour, we cannot 
expect an asymptotic expansion in a usual asymptotic scale for all these sequences. 
We provide only an asymptotic expansion in the mean, that is for the Cesaro mean 

1 ^ 

if we study the sequence (««)■ The way we will follow is based on elementary 
linear algebra. This is natural because radix-rational sequences are defined by 
linear recurrences and after all the mean asymptotic behaviour of a radix-rational 
sequence originates in the asymptotic expansion of a classical rational sequence, as 
we shall prove. Eventually we obtain a general theorem, valid for all radix-rational 
sequences. Evidently the asymptotic expansions which result form this theorem are 
not significantly different from those of our predecessors. Roughly speaking, they 
have the form 

J2 N" logi(iV) ^a;LiogB ^i^S^,i,^{logB N) + 0{N^'), 

a>a,,l>0 u 

where the w's are some complex numbers of modulus 1 and the ^'s are 1-periodic 
functions. The true difference lies in the fact that the theorem has a precise frame- 
work. 

The organization of the paper is as follows. In Section[2]we present the definition 
of a radix-rational series and a linear representation L, {Ar)o<r<B, C of such a 
series. We expose the link with the idea of a rational formal power series in non- 
commutative variables from the formal language theory. Such a series is the main 
object of study and a radix-rational series is only a metamorphosis of a formal 
power series through the interpretation of integers as their radix expansions in a 
given base. We address the problem of the asymptotic behaviour for K large of a 
vector- valued running sum 



Sa'(x) = Y. ^-^' 



\w\=K 

which depends on the length K of words and on a real number a; G [0,1] through the 
interpretation of words as B-a.ry expansions of real numbers from [0, 1). We give 
technical notations and hypotheses useful for this study. Particularly we recall the 
idea of joint spectral radius for a family of matrices, because the joint spectral ra- 
dius A* of a linear representation governs the precision of the asymptotic expansion 
we will obtain. Next we divide the running sum Sk^x) by its dominant behaviour. 
This gives a vector- valued function F/f (a;) and in the limit a functional equation 
appears. Surprisingly we have a dilation equation and such equations have been 
deeply studied because they are useful in the theory of wavelets and in the theory 
of refinement schemes. 

Section[3]is the basic part of the paper. TheoremlXl shows that the sequence (Fk) 
converges uniformly towards the unique solution F of the dilation equation. More- 
over in the special case where C is an eigenvector of a certain matrix Q relative to 
an eigenvalue with modulus p the speed of convergence is essentially 0{{X^/p)^). 
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In Section [4] we elaborate on this point by considering generalized eigenvectors. 
The main result is Theorem [BJ which gives an asymptotic expansion for Sk{x) 
from the linear representation. The technical hypotheses used as we work out the 
theorem reflect in the error term of the expansion which is essentially 0{X^). This 
asymptotic expansion has coefficients which are solutions of some dilation equations. 

Section [5] translates Theorem iBl for rational formal power series into Theorem [Cl 
for radix-rational sequences. So we obtain an asymptotic expansion for the running 
sums of a given radix-rational sequence (or for its Cesaro mean, it is the same). The 
solutions of the dilation equations become functions which are practically periodic 
with respect to the logarithm of the index. The theorem even though it is general 
may give an obvious formula. A typical example is the Thue-Morse sequence for 
which we conclude that it is 0(1). This flaw is quite normal because the running 
sum of the Thue-Morse sequence is almost the sequence itself. We illustrate the 
result with the study of the periodic function which appears in the asymptotic 
expansion for the Rudin-Shapiro sequence and specifically of its symmetries. We 
compare our work and the work of Dumont where a similar result appears. At the 
end we emphasize the fact that in full generality the coefficients of the asymptotic 
expansion are not periodic functions but only pseudo-periodic functions. 

To summarize, this paper gives a theorem, which is new and general, about the 
asymptotic behaviour of radix-rational sequences. Moreover it shows a link with 
the well developed domain of dilation equations. 

2. Radix-rational series and dilation equations 

2.1. Radix-rational series. For an integer n whose binary expansion is the word 
riK-i ■ ■ ■ no^ with figures nn-i^ . . ., ng in the set B = {0, 1}, the number of ones is 



k>0 



and the generating function is 



52(^) = 5]s2(ri)z". 

n>0 

It is easy to verify the following relationships 

VneN>o, S2{2n) ^ S2{n), S2i2n + I) = 82(71) + I. 

The constant sequence e = 1 satisfies obviously 

Vn G N>o, e(2n)=e(n), e(2n + 1) = e(n). 

It appears that the sequences S2 and e generate a Z>o-module which is left stable 
by the following linear right action of the monoid of words B*: the image of a 
sequence u under the action of a figure r e {0, 1} is the sequence u.r = {u{2n + r)). 
The example leads to the following definition. We consider a radix B > 2 and 
the associated alphabet B — {0, . . . ,B — 1}. For a semi-ring IK, the monoid B*, 
equipped with the concatenation, operates on the K-module of all sequences with 
values in IK by the map {u{n))„>o '~^ {u{Bn + r))„>o for each figure r in B. 

Definition 1. A generating Junction U{z) — X]n>o^"-^" whose coefficients take 
their values in a semi-ring IK is B -rational (or B -recognizable, or B -regular) if 
the ^.-module generated by the sequence under the right action of the monoid of 
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words B* is of finite type. It is radix-rational (or radix-recognizable, or radix- 
regular) if it is B -rational for some radix B. The sequence of coefficients (m„) is 
termed in the same manner. 



Allouche and Shallid ( 19921 ) named such sequences of coefficients "fc-regular se- 



quences" (fc is the radix). The adjective "regular" is reminiscent of "regular lan- 
guage" for a computer scientist but it is meaningless for a mathematician. The 
expression "recognizable sequence" has the same flaw, even if the preceding defi- 
nition would be more satisfying with recognizable in place of rational for a com- 
puter scientist. To the contrary "radix-rational sequence" make sense for both 
computer scientist and mathematician and we adopt this terminology. The name 
radix- rational series has the merit to remind that these series are a generalization of 
classical rational functions. These one admit a representation with only one square 
matrix, while _B-rational series admits a representation with B square matrices, as 
we shall see just below. 

Radix-recognizable series are a metamorphosis of recognizable formal series. In- 
deed, let 6""^ , . . ., e^ be a family of sequences which generates a module left stable un- 
der the action of B* and containing the sequence of coefficients u of a radix-rational 
series U{z). Such a family will be called a generating family for the sequence u in 
the sequel. Each of the sequences eKr, 0<r<B,l<j<d, writes as a linear 
combination e^.r = ^^ Or^ije^, not necessarily in a unique manner. The sequence u 
writes u = ^^ Cie', again not necessarily in a unique manner. 

Definition 2. The row vector L — (e^(0)), the matrices Ar — (ar.ij), < r < B, 
and the column vector C = (q), all together, is a linear representation of the 
series U{z) or of the sequence u. 

Besides the linear representation defines a recognizable series over B* with coeffi- 
cients in K in the sense of formal language theory. A formal series, usually written 

S = ^ {S,w)w, 
is an a pplication from B* into K, which associates to a word w the coefficient 



{S,w) (jBerstel and Reutenaueij . Il988t ISakarovitchl 120051 ). Here we consider the 



formal series S defined by 

[o,w) = -L^tui • ■ • A^j^ C 

for a word w = Wi . . .wk- The sequence u is obtained by restricting the formal 
series to the radix B expansions of integers. 

2.2. Notations and hypotheses. The data is a linear representation L, A ^ 
{Ax)x£Xi C of dimension d with a finite alphabet X of cardinality greater or equal 
to 2. The coefficients of the matrices are taken from the field of complex num- 
bers in the most general case. Up to a permutation, X may be taken equal to 
{0, 1, ... ,i? — 1} for some integer B > 2. The linear representation defines first 
a rational formal series S and second a radix-rational series U, by restriction to 
the radix B expansions of the integers. To abbreviate, we write A^j = A^^ • • • A^^^ 
for w = wi . . .WK, hence the writing {S,w) = LA^C. We denote by u the se- 
quence of coefficients of U , which means that we have u{n) = {S, w) for the integer 
n = {uK-i ■ ■ ■ «o)s whose radix B expansion is the word w = uk-i ■ ■ ■ no. 
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For such a given C-rational formal series S we consider the running sum over 
words w of length K submitted to the condition that the real number with radix B 
expansion (O.w)b is not greater than x in [0, 1] 

\w\=K 

for all nonnegative integer K and we want to estimate its asymptotic behaviour 
when K tends towards +00. With Q — Aq + Ai + ■ ■ ■ + Ab-i, it writes 

(1) Sk{x)= Y, Ar,Q''-^C+ Y A^^Ar^Q^'-^C + ^ A^.A^.Ar.Q'^-'^C 

ri<xx r2<X2 r3<X3 

+ ■ ■ ■ + / ^ AxiAx2 ■ • • Arj^O 

Tk<Xk 

if X admits the radix B expansion x = (O.X1X2 • • • )s- The formula renders evident 
the following lemma. 

Lemma 1. With Q = Aq + ^1 + ■ ■ ■ + Ab-i, the sequence of running sums (Sk) 
satisfies the recursion 

Sk+i{x)= J2 Ar.Q^'C + Ax.SKiBx-xi), 

ri<xi 

where xi is the first figure in the radix B expansion of x in [0, 1), with So(a;) = C. 

The matrix Q is the essential component which governs the mean asymptotic 
behaviour of S first, and of the sequence u next. Lemma [T] leads us to consider the 
powers Q^ of the matrix Q. More precisely, we look at the dominant term (from 
the asymptotic point of view) in the vector Q^C and we will use the following 



hypothesis (named (adt) for asymptotic dominant term). 

Hypothesis (adt) The sequence (Q^C) admits the asymptotic expansion 

(2) Q^C = R{K)V + 0{R'{K)) 

with 

and R'{K)^ RiK)o(^ 

Necessarily the vector V is an eigenvector of Q for the eigenvalue puo. 

In Section 21 we will expand the vector C over a Jordan basis, and this leads 



us to consider the following hypothesis (named (gev) for generalized eigenvector) 



This hypothesis leads to a more particular but more precise expression of Q^C 



than Hypothesis (adt) does. 

Hypothesis (gev) The family iV^^^)o<j<u is fuU rank and satisfies QV^"' — 

pwF(°) and QV^^^ = pujV^^^ + V'^^-^'> for j > 0, with p > and \uj\ = 1. 

At occasion, we will say that v is the height of the generalized eigenvector V^'^~^\ 
Formula (H]) use products of square matrices Ar, < r < B. A possible property 



to bound these products is the following (named (rsr) for rough spectral radius) 
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Hypothesis (rsr) There exists an induced norm || ||, and a constant A with < 
X < p such that all matrices Ar, < r < B, satisfy \\Ar\\ < A. 

We will tacitly use a norm on C'' and the induced norms on square matrices in order 
to guarantee the previous hypothesis. As a consequence we have HA^^ • • -^rjcH < 
X^ for every integer K. 



Hypothesis (rsr) will prove to be useful, but it is not sufficiently well designed. 
So we refine it as follows. We consider all the products A-^j for words w of a given 
length T and their norms. With the notation 

At = max P^,||^/^, 

■u) — T 



the joint spectral ra dius of the set ^r, < r < B, is the number (jRota and Strang! . 
19601 : lBlondell2008h 



A* = lim At. 

It is known that the joint spectral radius is not greater than any of the numbers At- 
Moreover A* is independent of the used induced norm. The following hypothesis 



(named (jsr) for joint spectral radius) is made to replace Hypothesis (rsr) 

Hypothesis (jsr) The joint spectral radius A* of the family of matrices {Ar)o<r<B 
is smaller than the number p, that is A* < p. 



It is known that A* is difficult to compute ( Tsitsiklis and Blondej 119971 ). For our 



purpose it is sufficient to find an induced norm and a T such that At < p- To the 
sake of clarity we will use a superior index to show the used norm if necessary, like 
Xj, if we use the absolute maximum column norm induced by the norm of index 1 
of C''. In the sequel, we will say that A* is attained if there exists an induced norm 
and an integer T such that A* = At. 

Besides, the joint spectral radius depends on the linear representation. To each 
representation is associated a finite dimensional vector subspace of the space of 
formal series left stable by the operators S i-^ S.r^^ = T^ ,(S,wr). The represen- 



tatio n is reduced if the subspace is as small as possible (jBerstel and Reutenaueii . 



19881 ). Evidently the smaller is the subspace the smaller is the joint spectral radius 



and all reduced representations provide the same A* because they are isomorphic. 
Hence it is better to always use a reduced representation and the joint spectral 
radius associated to a reduced representation depends only on the formal series; it 
is intrinsic. But it may be easier to compute the joint spectral radius for a non 
reduced representation, at the risk of obtaining a too large value. 

Example 1 [Rudin-Shapiro sequence] . The Rudin-Shapiro sequence may be defined as u{n) = 
(— 1)'^2 ;ii(") where e2;ii(n) is the numb er of (possibly overlapping;) o ccurrences of the pattern 11 
in the binary exp ansion o f the i ntege r n dBrillhart and Carlita . Il970f) . This sequence was defined 
independently bv lShapird 1 1195 j) and lRudinI 1119591) to solve a problem of optimality about the L°° 
norm of trigonometric polynomials with coefficients ±1. It is 2-rational: it admits the generating 
family {u{n), u{2n + 1)) and the reduced linear representation 

-(■ '). ^..-{i ;). ^^-{i \). c-[i 

The matrix 
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has two eigenvalues it\/2 which have the same absolute value. The vector Q^C has expression 



iK/ 



Q'^C: 



2K/2 ^ (2 + 2i/2) + (_i)A-(2_2i/2) 
21/2(1 + (-1)^') 



and H ypothesis ] (apt) | is not satisfied. The sequence is 4-rational too with representation IIAllouche and Shallij . 
l2003r) . relative to the generating family («(n), «(4n + 2)), 



(^ o, ^^-^-(l ; ' 4^^-^oA.^ ; V 



^"'-^-^"-^; ; . A?^Ai^(\ M, C^(\ 



The matrix Q^^' = A^ + AqA^ + AiAq + A\ = {Aq + A-i)'^ = Q^ has a dominant eigenvalue p = 2 
and it is evident that Hypothesis |(adt)| is satisfied for this linear representation. Because Aq 
and Ai maps each vector of the canonical basis onto a vector of the canonical basis or its negative, 
we have \t = 1 for every T > 1 if we use a norm which gives the same value for both vectors of 
the canonical basis. Hence Hypothesis | (rsr)] is satisfied for the radix 4 representation above and 
Hypothesis | ( jsr)| too because A* = 1. 

2.3. Self-similarity. For any positive integer K, let us introduce tiie function F^ 
from [0, 1] into C* defined by 

Lemma [1] translates into the equation 

We may consider the operator Hk of the space of continuous functions from [0, 1] 
into C* defined by 

^ ' r\<x\ ' 

(In all the paper x\ is the first digit of x G [0, 1).) Equation Q rewrites 

Fx+i(a;) ^ LkFk{x). 



According to Hypothesis (adt) the sequence of operators Lk converges weakly 



towards the operator C defined by 

(5) L^(x) = — y Ar^ ^-—A^,^{Bx-x{) 

r\<x\ 

and we will first study the equation £$ = $. 

In the sequel, the following system of equations, whose unknown is a function 
from the segment [0, 1] into the space C*, 

- $(0) = 0, $(1) = F, 

- for every figure r of the radix B system and for x in [r/_B, (r + 1)/_B), 

(6) $(a;) = — VA.,V^+ — Ar^(Bx-r) 

puu ^-^ puj 

ri<r 

will be named the basic dilation equation. 

Proposition 1. Let L, {Ar)Q<r<B, C he a linear representation of dimension d 



for the radix B . Under Hypotheses (adt) and (rsr) the problem 

- ^ is a continuous function from the segment [0, 1] into the space C^, 
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- ^ is a solution of the basic dilation equation, 
has a unique solution F. 

Before we go further we have to do a simple remark. If we multiply the square 
matrices A^, < r < B, hy a nonzero scalar a, the eigenvalue puj is multiplied 
by a, the running sum S/f (x) and the dominant asymptotic behaviour R{K) are 
multiplied by a^ . It follows that F^ and the operator Ck are unchanged. In the 
same manner the operator C is not modified and the solution of the fixed point 
problem described in the previous proposition remains the same. This permits 
us to consider that the modulus 1 number tu is equal to 1 in order to prove the 
assertions. Even we may assume that puj is equal to 1 but we prefer to keep in 
mind the modulus p of the eigenvalue puj. (Nevertheless see the next subsection.) 

Proof. According to the previous remark, we may assume a; = 1 by multiplying all 
the square matrices A^, < r < B, hy the conjugate number uj. 

The space of continuous functions from [0, 1] into the space C^, equipped with 
the norm of the maximum H^Hq^ — max^^ ||<I>(a;)||, is a complete normed space. 



(Recall that Hypothesis (rsr) assumes that we have chosen a norm on C''.) The 
continuous functions $ which satisfy $(0) = and ^{l) — V are the elements of 
a closed, hence complete, subspace C of this complete space. The equation of the 
problem appears as a fixed point equation $ = £$. It is sufficient to see that the 
subspace C is left invariant by C and that £ is a contraction to prove the assertion. 
We must verify that for a continuous $ in C the function ^ — £$ is a member 
of C. According to the piecewise definition of 'J, we have to consider the left and 
right limits of ^ at the points r/B for < r < B. The definition of ^ and the 
continuity of $ give immediately 

*(0) = -Ao<^iO) = 0, 
P 

*(i-o) = - Y. AsV + -Ab'V = -qv = v, *(i)-y. 



P oil^B' P P 



and for < r < B 



^ Q<s<r 

B p ^-^ p p ^-^ 

^ 0<s<r-l ^ ^ 0<s<r 

The constraints are satisfied and the subspace C is stable. 

If we have two functions $1 and $2 in C, let ^1 and ^2 be their images by L. 
From \I'i(x) - ^2(2^) = (l/p)^r(*i(-Ba; - r) - '^i{Bx - r)) follows the inequality 
11^1 — ^2 II < (-^/p) 1 1*5*1 ~ *&2|| because the norm of each matrix A^ is bounded 



by A. As a consequence of Hypothesis (rsr) the operator £ is a contraction. D 



In fact Hypotheses (adt) and (rsr) are not necessary to conclude that the 



solution of the basic dilation equation is unique. A drop of regularity is sufficient. 

Lemma 2. The basic dilation equation may have only one solution under the sole 
hypothesis p > 0, |ci;| = 1 if we ask for a function continuous on the right. 



10 
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Proof. For a _B-adic number, that is a number x whose radix B expansion is finite, 
the dilation equation gives the value F{x) because the recursion formula provided 
by the dilation equation terminates on the basic case x = 0. But the set of _B-adic 
numbers is dense in [0, 1] and the function is assumed to be right continuous at 
every point. Hence it is completely determined. D 

Example 2 [Regular self-similar functions]. Often the functional equations like ^ are pre- 
sented as leading necessarily to chaotic functions, but some very regular functions may satisfy such 
a self-similarity property. As an example consider the sequence with takes value 1 on all the multi- 
ple of 3 and for the others integers. It is rational in the classical sense and its generating function 
is rat ional with poles which are all roots of the unity. As a consequence (AUouche and Shallig, 
I2OO3I . Th. 16.4.3, p. 446) the sequence i s rational with respect to every radix. Here is a binary 



linear 


representation dSakarovitchl. |2005|. d. 


6), 






/ 1 

L=(l 0), Ao= 

V 1 




1 



\ / 1 \ 

, Ai = 10 
/ V 1 / 


/ 1 \ 
C= . 
V / 


The system satisfied by F is 










Fi(x) = -Fi(2x), 




' Fi(x) = i + lF2(2x-l), 




■ 


F2(x) = -F3{2x), iorO<x<-; 
2 2 


< 


^2^ = ^ + ^^1(2^-1). 


for < x < 1. 




^ F3(x) = iF2(x), 




^3W = i + iF3(2x-l). 





W e find immediately the solution F(a:) = ( x x x ) . 

Di nsenbacher and HardinI ( 19991 ) have studied the distributional solutions with 
a bounded support of a dilation equation. The idea is to use an antiderivative of 
sufficiently high order n. This makes contracting the operator behind the equation. 

2.4. Wavelets and refinement schemes. The basic dilation equation enters into 
the domain of what is called a two-scale difference equation, namely 



V{x) = ^ Cn(p{2x - n). 



(For the sake of simplicity, we limit ourselves to radix _B = 2 in this subsection.) 
These equations have been heavily studied because they appear in the theory of 
wavelets to define a scale function and in the theo ry of refineme nt schemes of com - 
puter graphics to define a refinement function. iDaub echies a nd Lagarias ( 1991) 
provi de an expository of their occurrences and a bibliography. See also (JHeil and Colella , 
19961) . 

For compactly supported wavelets the previous sum is finite and the equa- 
tion may be rewritten in a way wh ich looks lik e our dilation equa tion. We fol- 
low (JDaubechies and Lagariasl Il992l . p. 1036) or (JDaubechiesJ . Il992l p. 235). The 
equation 



ip{x) 



N 

E 



Cn(p{2x — n) 



is translated into an equation 

ro$(2.T) if 0<a; < 1/2, 

Ti<^{2x - 1) if 1/2 < a; < 1, 

with a vector-valued unknown function 



$(a;) 



<^{x) = {lp{x) (^(x-fl) (p(a; + 2) ... Lp{x + N - I) ^ 
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and square matrices Tq — {c2i-j-i)i<i,j<n , Ti ~ [c2i-j)i<lj<n (it is assumed c„ = 
for n outside [0, N]). Scalars c„ are constrained by X^m C2n — X^n C2n+i = 1 and 
the number 1 is an eigenvalue for both Tg and Ti, and for M = {c2i~j)i<i,j<N-i- 
We consider a right eigenvector W = {wi)i<i<N-i for Af (an additional condition 
imposes that 1 is a simple eigenvalue, so there is essentially only one possibility 
for W). We extend W hy wq = 0, wn = and define V — {vi)i<i<N by Vi = Wi-i 
for 1 < « < TV, so that W is an eigenvector for Tq and V is an eigenvector for Ti 
relative to the eigenvalue 1. Boundary conditions are added to the equation, namely 
$(0) = W, $(1) = V, which guarantee the well definition of the equation and the 
continuity of the solution ip (under some conditions on the c„ 's) . 
Besides our dilation equation writes 

ToF{2x) if < a; < 1/2, 



^^^' \ ToV + TiF{2x-l) ifl/2<j;<l, 
where matrices Tq and Ti are defined by 

To = — Ao, Ti = — Ai, 

pu) puj 

and V is an eigenvector for To + Ti relative to the eigenvalue 1. Here the boundary 
conditions are F(0) = and F(l) = V. 

Evidently both equations are very akin. The two-scale difference equation for 
the scale function of wavelets is homogeneous while our dilation equation is not, 
but thei r linear parts are t he sam e. Inhomogeneous dilation equations have been 
studied ( Strang and Zhoul . Il998r ). because they are useful in the construction of 



wavelets on a finite interval and of multiwavelets. Nevertheless this contrast is an 
illusion: if we extend the function F as a continuous function over the whole real 
line by making it constant on (— oo,0] with value and constant on [l,-|-cx)) with 
value V, the dilation equation rewrites 

F(a;) ^ ToF(2.x) + TiF{2x - 1) 

for X real. The writing of the dilation equation by cases is more concrete, but the 
homogeneous version above is more compact and more practical for proofs. So we 



will made use of the following convention (named (hec) for homogeneous equation 
convention) at occasion. 

Convention (hec) The solution of a basic dilation equation is extended to the 
whole real line as a continuous function constant on the left of and on the right 
of 1. 

The eigenvalue 1 appears in both cases. The boundary conditions are not of the 
same form and the matrices To and Ti for the wavelets have a very special struc- 
ture, while the matrices Aq and Ai of a linear representation are not constrained 
in our study. As a result, even if the computations are not exactly the same, the 
ideas which work for wavelets work too for ratio nal series. For example, the basic 
idea of t he cascade algori t hm (ipaubechiea . ll992l . § 6.5, p. 206) or of the refinement 



schemes (|Dvn and Levinl . 120021 ). which computes the value of a scale function for 
dyadic numbers, applies here. We have yet used it in Lemma [H and evidently to 
draw the pictures in the paper. In the same manner, the key point for the existence 
and uniqueness of the solution of these two-scale difference equ ations is the occur- 
rence of a contracting operato r (iDaubechies and Lagariaslll99ll . Sec. 4). The same 



idea have appeared before in (JHutchinsoiJ . Il98ll ) which describes a construction of 
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self-similar parameterized curves (and the construction of the sequence (Gk) in 
the proof of Lemma [3] below is of the type described in its § 3.5). 

3. Basic limit theorem 

3.1. Uniform convergence. We are now in position to prove a first result about 
the asymptotic behaviour of the running sums Sk- 



Proposition 2. Under Hypotheses (adt) and (rsr) , the sequence {Fk) converges 
uniformly towards the function F defined in Proposition d 

Proof. To obtain a uniform convergence, we take care that all the big oh with 
respect to K are uniform with respect to x. Also we may assume w = 1 as in the 
proof of Proposition [TJ 

We first note that the sequence {Fk{x)) is uniformly bounded. Actually let e > 
chosen such that 

p I — e 
There exists a Kq such that for K > Kq we have the inequality 



1R{K+1) _^ 



< e. 



P R{K) 

The triangular inequality applied to the right member of Formula Q provides an 
inequality 

||FK+ilL<fc(7+||FK|L), 
where 7 is a constant. By induction we obtain 

I|FkIL<Y^ + I|FkoIIoo 

for K > Kq. Hence the sequence (||F7<||j^) is bounded. 
Formula (|4]) and Formula ([6]) (applied to F{x)) provide 



(7) FW.)-FW.i:A..(^^e'-C-V) 



ri<xi 



A ( ^(^) , Fk(Bx - xi) - -F(Bx - xi 



The writing Q^C = R{K)V + R{K)0{1/K) gives 

Fk+i{x)-F{x)^ J2 Ar,o(^)+A,,FK{Bx-xi)0 



+ -A^, {Fk{Bx - xi) - F{Bx - xi)) 
P 

and next, with the uniformly bounded character of {Fk{x)), 

||FK+i-F|U<0(^-i)+^||F,,-F|U. 

By induction, we find \\Fk — F||^ = O (l/K). If the constant implied in the big 
oh of the last but one formula is 7, we may take 27/(1 — X/p) for the constant in 
the big oh of the last formula. D 
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We may be more precise about the speed of convergence, but we content ourselves 
with a particular case (which will prove to be basic). 



Corollary 1. Under Hypotheses (adt) , (rsr) if C is an eigenvector associated 
to the eigenvalue pus, the speed of convergence is \\Fk ~ F||j^ = 0{{X/ p)^). 

Proof. We do not assume here \uj\ = 1 hence p must be replaced by pto. By 
substituting Q^C = R{K)V + 0{R'{K)) into Formula dT]), we obtain 



(8) 



' A'+l 



<o( ^(^) 

- \R{K + l) 



1 

puj 



O 



( R'{K) 



X 

p 



K 



But if C is an eigenvector, we have R{K 



\R{K + l) 
l)/R{K) = puj and R'{K) 



Formula ^ becomes ||F 



K+l 



< (A/p) \\Fk 



and 
hence the conclusion. 

D 



3.2. Basic theorem. Let us assume that we group the letters into pairs, which 
means that we consider only words of even length. We obtain a new formal series 82- 
It is rational and admits a linear representation whose square matrices are the B^ 
products ArAs with Q < r,s < B. The associated matrix Q becomes Q2 = Q^ ■ (See 
Ex. [TJ) The sequence of running sums (S^) is changed into its subsequence {S2K) 
and it is the same thing for the sequence {Fk)- In the same manner we may group 
the letters by T for a given T. This leads us to consider the subsequence [Ykt) 
and the power Q^ . 

In order that Proposition [5] applies to a power of matrix Q, we have to consider 
all products Ayj for words w of a given length and their norms. Hypothesis (rsr) 
will be satisfied for some power Q^ if we impose Hypothesis (jsr) 
chosen such that A, (1 + e) < p and T is sufficiently large, we thus have \\A. 
(Xt)^ < Xj{l+e)'^ < p^ for all words w of length T. Hence the subsequence {Ykt) 
is convergent. We will show that not only this subsequence is convergent but the 
sequence (F^) is convergent. 



If £ > is 
< 



Theorem A. Under Hypotheses (adt) and (jsr) , the sequence (F^) converges 
uniformly towards the unique continuous solution F of the basic dilation equa- 
tion ([5]). 

Proof. As we have explained just before the wording of the theorem, the sequence 
{Fkt) converges uniformly to a function F" for some T. Let us consider the 
sequence (F^T+i)- According to Lemma [U wc have 



R{KT + 1 
hence 



R{KT+1) 



R{KT+1) 



' KT+l 



(x) 



R{KT 



1) ^ 

ri<x-i 



Ar^ {R{KT)V 



0{R'{KT))) 
R{KT) 



' KT 



{Bx — xi). 



R{KT+1) 

Because (Fkt) converges uniformly towards F°, we see that (Fkt+i) converges 
uniformly towards F^ = /CF". Repeating the argument, we conclude that each 
subsequence (Fkt+s) converges uniformly towards F" = £^F° for s > 1. But each 
of these functions satisfies F* = £-^F'', because F° is a fixed point of C'^ . Since 
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the unique solution of this equation is F'', aU these functions are equal. Henceforth 
the sequence (F^) is uniformly convergent towards F". Because of the equality 
po _ pi ^j^g funct^JQn po jg ^^^ unique solution of the equation $ = £$. (The 
uniqueness of the solution for the last equation is guaranteed by Lemma [2]) D 

The following case will prove to be useful in Section[H where it will be extended. 



Corollary 2. Under Hypotheses (adt) and (jsr) , ifCis an eigenvector associated 
to the eigenvalue puj, the speed of convergence is \\Fk — F\\^ = 0{{X/p)^) for every 
A > A*. If X* is attained we may replace A by X^. 

Proof. Corollary [1] gives the speed of convergence for the subsequence (Fkt) and 
we obtain \\Fkt — FIIoo — 0[rff) with r^ = XTj./ p^ because we use the power Q^ 
in place of Q. This gives ||Fkt+s — F||^ — 0{r^) for < s < T because the 
operator C is continuous and because there is a finite number of subsequences to 
consider. We obtain \\Fk — F||^ = 0{rj,' ) and the formula above because T was 
chosen such that At < A*(l + e) for an arbitrary e > 0. In case A* is attained the 
positive number e is useless. D 

Example 3 [Worst mergesort sequence]. Mergesort is a comparison based sorting algorithm 
which uses the recursive divide-and-conquer strategy. The Ust to be sorted is spht into two Usts 
of almost equal size; both are sorted by mergesort (there is nothing to do for a list with only one 
item); both sorted lists are merged. Taking into account the number of comparisons, the cost of 
mergesort for a list of n items is 

wher e m^ is the cost of merging two sorted lists with [n/2j and [n/S] elements IIFlaiolet and Golinl . 
119941) . The cost of mergesort in the worst case (that is m„ = n — 1) writes ii„ = X]fc=i ri°g2 ''^l 
and the sequence (fn) defined by no = and d„ = [log2 n] for n > 1 is a 2-rational sequence 
which admits the representation, relative to the basis (i'(n), v(2n), v(2n + 1), ti(4n + 2)), 



Ao 






V 1 ) 

The matrix Q = Aq + Ai has two eigenvalues, namely 2 and 1, which are double. A computation, 
based on the Jordan reduced form of Q, gives 



( ° 


-1 


-1 


-1 \ 




/ 


1 


1 




2 



1 
1 


1 




, Ai = 




1 


-1 
-1 


V 








1 / 




V 1 


2 



1 \ 




( ^ ^ 


-1 

-2 


, c = 






3 / 




V ; 



Q'^C-- 



2^K 



/ \ 




/ \ 




/ -1 \ 




/ 1 \ 





+ 2^' 





+ K 


1 


+ 





-1 




2 




1 




-2 


V 1 / 




V -1 / 




V -1 / 




V 1 / 



and Hypothesis | ( apt)] works with 

p = 2, u] = l, R{K) = 2'^K, 



y = ( -1 1 )' 



R'{K) 



With the absolute maximum column norm, we find that the maximal norm of the matrices Ay^ 



S.W 



with w of length 4 is 9. This gives \)^' ~ 1.73 and Hypothesis | (.isr)| is satisfied. As a consequence 
(Fk) converges uniformly towards F which is nothing but 

F(a:) = ( -x x )'" . 

3.3. Holder property. Let us recall that a function / from the real line into a 
normed space is Holder with exponent a if it satisfies ||/(j/) — /(a^)|| ^ c\y — a;|" 
for some constant c. Under the hypotheses of Theorem [X] the function F is not 
only continuous but Holder. Such a result is classical in the study of dilation 
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equat ions and the proof below is a replay of point (7) in the proof of Th. 2.2 



from (JDaubechies and Lagariaa Il992l ) . See also (|Riou]| . Il99 



i 



Lemma 3. Under Hypotheses (adt) and (rsr) the function F is Holder with 
exponent logg{p/\). 

The inequality A < p is assumed, but we haye also p < \\Q\\ < BX by triangular 
inequality. This giyes bounds < log^{p/X) < 1 for the exponent, as expected. 

Proof. Let us introduce the sequence G^ defined by Go (a;) = xV and Gk+i = 
CGk, where C is the operator defined by Eq. ([5|) and used in the basic dilation 
equation © . The function Gq is linear and thanks to the properties of the opera- 
tor £, all the functions Gk are piecewise linear and continuous. More precisely Gk 
is Hnear in each of the interval [k/B^, {k + l)/B^], 0<k< B^ . 

Let s and t be two real numbers which are situated in one of these intervals 
[k/B^, {k + l)/B^), 0<k<B^, and satisfy 1/B^+^ < \t - s\ < l/B^. Their 
radix B expansions have the same first K figures, but the {K + l)th figures are 
different. We may write s — (O.ui . . . ukSi . . ■)b, t — (O.ui . . . uxti . . ■)b and s' — 
(O.Si . . .)b, t' — (O.ti . . .)b- According to the definition of £, wc have 

1 





Gi, 


■it)- 


-Gk 


(s) 


{puj) 


K ^"1 ■ ■ 


Auk 


it'-.. 


/)y. 


Consequently we 


obtain 
















\\Gk 


w- 


-GKis)\\ 


< 


'^V \ 


V\\ < B 


\\v\\ 


fBX^ 
\ P y 


K 

) 1' 



using the hypothesis 1/_B^+^ < |i — s|. Because Gk is piecewise linear this in- 
equality is valid for all pairs {s,t) taken from [0,1]. Let us insist on that point. 
The function A : {s,t) ^^ \\GK{t) — Ga'(s)|| /\t — s\ is defined on the square [0, 1]^ 
minus its diagonal. It gives the mean speed of the parameterized curve Gk on 
the interval whose ends are s and t. According to the previous computation the 
quantity B \\V\\ {BX/ p)^ is an upper bound for the mean speed A on each square 
[k/B^ , {k + 1)/B^) (minus their diagonals), < fc < B^ , because Gk is piece- 
wise linear and A is constant on such a square. Again because Gk is piecewise 
linear these particular squares give the larger value of A. Henceforth the previous 
upper bound is valid on the entire square [0, 1]^ (minus its diagonal). 

Let X and y be two real numbers from [0, 1] which satisfy 1/B'^'^^ < \y — x\ < 
l/S" for some integer v. The sequence (Gk) converges uniformly towards F and we 
may guarantee ||Gi, — F||^ < "f[X/ p)" for some constant 7 because the operator C 
is a contraction with ratio X/ p. Using the triangular inequality we obtain 

||F(j/)-F(x)||<27Q) +B\\V\\(^^^ \y-x\<{2^ + B\\V\\)[j^ 

the last inequality coming from the hypothesis \y — x\ < 1/ B'^ . But the inequality 
\/ B'^ < B\y — x\ provides an inequality of the desired form with the help of the 
formula X/p = [1/ Bf°^B{p/>-) , D 



Proposition 3. Under Hypotheses (adt) and (jsr) the function F is Holder with 



exponent log^(/9/A) for every A > A*. Moreover A = A* may be used if A* is 
attained. 
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Figure 1. Limit probability distributions which come from a 
BernouUi process described in Ex. |4l 



Proof. If A > A* there is a T such that A > At and we may assume A = At for 
the proof because the smaller is A, the stronger is the constraint imposed by the 
Holder property (the exponent logg{p/X) is a decreasing function of A). We replace 
the radix B by B^ . This changes p into p^ and A into A"^, but loggT{p'^ /X^) = 
\ogg{p/ X) remains the same and the previous lemma gives the conclusion. D 



Example 4 [Billingsley's distribution function] . iBillingsIevI l|l995l . Ex. 31.1, p. 407) studied 
the random variable X = X]n>0 ^i/2" where Xn is the result of a coin tossing with probabili- 
ties po and pi for Xn = and Xn = 1 respectively. This defines a rational series with dimension 1, 
radix 2 and a linear representation 



( 1 ) , Ao = { PO ): 



(Pi ). c={ 1 ) 



with < po,pi < 1, PO +Pi = 1. We have Q = { 1 ), p = l,V = { 1 ), R{K) = 1, R'{K) = 0, 
A, = max(po,pi). The distribution function F is the limit function of Theorem 1X1 and it is Holder 
with exponent l og9(l/A«). We i llustrate the example (Fig. [T] left-hand side) with po = 1/4, 
pi = 3/4 as in l lBillingslevl Il995l. p. 408) and the exponent is log2(4/3) ~ 0.415, and (Fig. [T] 
right-hand side) with pp = 1/5, pi = 4/5 and the exponent is log2 (4/5) ~ 0.322. Following the 
same way as in IjDumas et al.l . 120071 . § 4.3), it is possible to show that, assuming po < pi, at every 
dyadic point the best Holder exponent is log2(l/po) on the right-hand side and log2(l/At) on the 
left-hand side. Except in the case po = Pl = 1/2, this gives log2(l/po) > 1 and this explains the 
horizontal tangents on the right-hand side that we see on the pictures. But we will not elaborate 
upon this point because the argume nt is very simple and ass umes a linear representation with 
nonnegative coefficients while llDaubechies and Lagariasl Il992l) has provided a general approach 
to this subject. 



4. Asymptotic expansion 

To a linear representation L, {Ar)o<r<B, G, we associate a vector-valued se- 
quence of running sums 

Sk{x) - Y. ^-'^■ 

\w\=K 
{0.w)b<x 

We want to derive an asymptotic expansion for the sequence (S^) and as a by- 
product an asymptotic expansion for the sequence {Sk) = (LSk)- Theorem lAl and 
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its corollaries provide a one term expansion 

Sk{x) = R(K)F{x)+o{R{K)). 

In order to obtain a more precise expansion, we generalize the result of Corollary[2l 
which deals with eigenvectors of Q. We find a basis which reduces the matrix Q = 
'^o<r<B ^r to its Jordan normal form and we decompose the column vector C on 
this basis. The sequence (Sk) appears as a sum of sequences associated to the 
vectors of the basis. We have to discuss according to the eigenvalue of Q relative to 
each vector. If the modulus of an eigenvalue is larger than A* we have to consider a 
generalized eigenvector or a Jordan vector and this will be made in the next section. 
After that it remains only to consider the case where the modulus is less or equal 
to the joint spectral radius A*. This case produces a noise which will enter in the 
error term of the asymptotic expansion. 

4.1. Jordan vector. We use a family of linearly independent vectors (V^*--''')o<j"<iy 
which satisfies QV'^^^ = puV^^'^ + V^^'^^ for j > and QV^ = pujV'^^^ with 
p > A* and \uj\ = 1, that is we assume Hypothesis (GEv)| with Hypothesis (jsr)[ 
As a consequence Q induces on the vector space generated by V^^"^ , V''-^^ , . . ., Ia^^^ 
the usual Jordan block of size v, 



(9) 



'J puj 



( p^ 


1 

pu 


1 


\ 


K 






1 

puj ) 



This gives immediately 

(10) Q^i/('''i) = (^_ \ (pc^)^-'^+ V(") + if^ {p^t-"^^y^^^ + • ■ • 

We claim that the running sum associated to the vector y(''^i) 

\w\=K 

admits an asymptotic expansion of the form 



Sk(x) 



K 

K-*+oo V J/ — 



^ypuf'^^'F^^Hx) + (^^2)(pc.)^-''+^F(^)(x) + 



K 



(puj) 



A'-lTp(i/-2 



Hx) + {pu;fF'^''-^\x)+erTOTK{x), 



where the F'^^'s are continuous functions from [0, 1] into C^ and the error term is 
O(A^) for every A > A*. 

The polynomial function K i~> (pcu)^^ Q^ V^'^^^' has a unique writing on the 
basis (, ), < j < i^, and we have necessarily F'^^'(l) = V^^^ for < j < i^. In 
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the same manner the uniqueness of asymptotic expansions shows that the family 
(F(-'))o<j<iy must satisfy the system of dilation equations, 

ri<xi 



(11) { 



F(''-2)(a;)+pcjF(''-i)(x) = J" Ar,V^''-^'> + A,,F^''-^'>{Bx - xi). 



ri<xi 

We obtain these formulae by substituting the asymptotic expansion into the func- 
tional equation of Lemma [T] Evidently the first equation of the system is nothing 
but the basic dilation equation ([6]). To obtain the proposition we have in mind, we 
will follow the same path as in Sections I3.1H3.3I For the sake of clarity we cut the 
proof into lemmas. 



Lemma 4. Under Hypotheses (jsr) and (gev) , System (jlip has a unique solu- 
tion (F^-'-')o<j<iy in the space no<j<,yC'^' where C^^' is the space of continuous func- 
tions from [0, 1] mto C^ which satisfy F(j)(0) = 0, F(j)(1) = V'-^'^ for 0<j<iy. 

Proof. It is possible to consider one equation at a time, but it is more enlightening 
to deal globally with the system. Let us introduce the matrix F(x) of type d x [0, ly) 
whose columns are the column vectors F^''^(x), F'-^^x), . . ., F'-'^~^\x). The system 
writes 

(12) ¥{x)Jp^ ^Y,ArV + A,,¥{Bx - xi), 

r<xi 

where V is the matrix whose columns are the column vectors V^'^\ V^^^', . . ., V^'^^^\ 
and Jp^ is the Jordan block defined by Eq. ^ . Because of the positivity of p, the 
matrix Jp^ is invertible and the system turns out to be a fixed point equation 
F = C¥, namely 

¥{x) = Y, ArVJ-J + A,,¥{Bx - xi) J-J 

r<xi 

for some operator £. Let /z lie between A* and p, that is A* < /i < p. Because 
the spectral radius of Jz} is l/p, there exists an induced norm which provides 
1/p < ll'^pw^ll < 1/m- Next there exists an integer T which gives A, < At < M- We 
consider for a while that the radix is B^ . From \\Axi \\ ||^p^"'^|| < /^ x 1/p = 1, we 
see that the right-hand side of the previous equation turns out to be the image of F 
by a contracting operator. We conclude as in the proof of Proposition [1] that the 
equation F = C'^¥ has a unique continuous solution F^^^ from [0, 1] into the space 
of matrices Mdx[o,z.)(C) which satisfies the constraints F(°)(0) = 0, F(°)(1) = V. 
Next we prove that F^^^ is the unique solution of the problem, in the same way we 
have finished the proof of Theorem [X] D 
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Lemma 5. Let Ak be the function from [0, 1] into C* defined by 

Ak{x) ^ ^f_ J (pc.)^-+iF(")(x) + ^J[^ [pu.y<-^+^F^'\x) + ■■■ 

for each nonnegative integer K , where F^'^' , . . ., F'"^"^) are the functions whose exis- 



tence and uniqueness are asserted by Lemma^ Under Hypotheses (jSr) and (gev) 
the running sum Sk associated to V^'^^^' writes Sk{x) — Ak{x) + 0{X^) for ev- 
ery A > A* and the big oh is uniform with respect to x in [0, 1]. Moreover if X^ is 
attained we may replace A by X^. 

Proof. With errori<-(a;) = Sk(x) — Ak{x) it is not difficult to obtain, by substitu- 
tion, the recursion 

error j<-+i(a;) = A^-^ eiroi k{Bx — xi). 

More precisely, we use the recursion of Lemmall] System ()lip . and the fundamental 
formula of binomial coefficients. For A > A* we choose an integer T such that 
A* < At < A. We readily obtain ||errori^:T|loo = 0(A^^). Next the recursion shows 
that we have ||errori<-T-(-s|loo = 0(A^"^^*) for < s < T, hence the result. If A* is 
some At, it is useless to consider a A > A*. D 

As in Lemma m we show that the functions F'^-'^'s are somewhat regular. 



Lemma 6. Under Hypotheses (jsr) and (gev)[ each of the functions F'^', < 



j < v, defined in Lemmal^is Holder with exponent logg{p/X) for every X > A*. // 
A* is attained, they are Holder with exponent logglp/X^). 

Proof. We define a sequence of parameterized curves in the space of matrices of 
type d X [0, f). The first one is /q = and the others are defined by the recursion 
Fk+i — CFk, where C is the operator we used in the proof of Lemma IH namely 

rK+l{x) = J2 Ar-yJ-J + A^.PKiBx - Xi)J-J. 
r<xi 

Each Fk is piecewise linear. 

We choose an induced norm and a T such that A* < At < p, and we consider 
that we use the radix B^ . As in the proof of Proposition [3] we observe that the se- 
quence (Fkt) converges uniformly towards the solution F = (F(^')o<j<;y described 
in Lemma |4l This permits us to show that F is Holder with exponent logg{p/XT). 
The argument remains the same: the maximal speed of a piecewise linear function 
is the greatest uniform speed on each piece and to compute a uniform speed it is 
sufficient to test two distinct points of the interval under consideration. D 

We summarize the obtained result. 



Proposition 4. Under Hypotheses (jsr) and (gev) , System pT|) has a unique 



solution (F'-'')o<j<iy in the space of continuous functions from [0, 1] into C^^P'"^-* 
such that F(j)(0) = 0, F(-')(1) = V^'^ for < j < ly. Each of the functions F^-?), 
< j < i^, is Holder with exponent logg{p/X) for every A > A*. The running sum 
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S_R- associated to V^"^^^' admits the asymptotic expansion 

for every A > A* and the big oh is uniform, with respect to x in [0, 1]. Moreover if X^, 
is attained, the functions F^^' 's, < j < v, are Holder with exponent \ogg{p/X^,), 
and we may replace A by A» in the asymptotic expansion. 

Despite of its technical aspect, the obtained result is very simple: each term of 
the expansion of Q^F*-'^"^-' is shaken by a somewhat regular vector-valued function 
which goes from to the vector which appears at the same place in the expansion. 

4.2. Noise. So far we did not take into account the case p < A*. Some experiments 
lead us to expect a chaotic behaviour for the running sums and we content ourselves 
with exhibiting a natural bound. 

First we consider an eigenvector V oi Q associated to an eigenvalue puj, with 
p > and \uo\ = 1. Note that this time the modulus p may be 0. 

Lemma 7. Let V be an eigenvector for Q with eigenvalue of modulus p > which 
satisfies p < X^. The sequence of running sums associated to V satisfies W^kWoq — 
0{X^) for every A > A*. //A* is attained, we may replace 0{X^) by 0{X^) in 
case p < A* and by 0{KX^) in case p = A*. 

Proof. The functional equation satisfied by the running sum associated to the eigen- 
vector V is 

8^+1(0;)= Y. Ar.Q^'V + A^.SKiBx-x,), 



that is 



ri<xi 



Let us assume first that we have a A such that p < X and \\Ar\\ < X for < r < B 
and for some induced norm, that is let us assume in a certain sense the converse of 



Hypothesis (rsr) We obtain 

\\Sk+i{x)\\ < p"" Y ^\\V\\+X\\SK{Bx-Xi)\\, 



ri<xi 



hence 



00 ■ 



\\Sk+i\\^<p''BX\\V\\+X\\Sk\ 

With ||S/i-||j^ = X^uk, we have uk+i < (p/X)^ B \\V\\ + uk- In case p < A, we 
have Uk < u^ + B \\V\\ /(I — p/A) and we conclude ||S/f ||^ = 0{X^). In case p = X, 
we have uk < uq + KB \\V\\ and we conclude ||Sif ||^ = 0{KX^). 

Second we consider the general case, which divides into two sub-cases. Let 
us assume that the joint spectral radius is smaller that all At, for every induced 
norm. For any A > A*, we may find, for a given induced norm, a T such that 
A* < At < A. We apply the previous argument to the sequence (Skt)- This 
gives IJSktIIoo = 0{X^^) because we have p < A* < At- Next, Lemma [1] gives 
llSifT+slloo = 0(A^'"^+'*) for < s < T, because there is a finite number of s 
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Figure 2. The image of the parameterized curve Sk from Ex. [5] 
grows hnearly with K. Here we see S3, Se, and S12. 



O(A^), and we conclude 



between and T. In this way, we obtain ||Si<-||j^ 
l|SKlloo = 0(A^)foreveryA>A,. 

Let us assume that A* is some At for some induced norm. In case < p < A* = 
At the same reasoning as above, with A* in place of A, gives ||Si<-||^ = O(Af^). In 
case p = A* we find ||S_r-t||oo = 0{KX^^). As above and again because there is a 
finite number of integers between and T, we obtain ||S/i-||j^ = 0{KX^). D 

The next example shows that the upper bound we have found in case p = A* 
has the right order of growth. 



Example 5 [Triangular tiling]. For a real •&, we consider the rotation matrix 

cos 1? — sin 1? 
sin 1? cos 1? 



R^ 



In the following we use Vq = Ei the first vector of the canonical basis and V^ = R^Vq. The 
example is based on the linear representation 

L=(l 0), Ao = i?_^, Ai=Ra, C = Vo. 

Because we use orthogonal matrices, the joint spectral radius is A* =1. Wc specialize the case to 
1? = 7r/3, so that we have p = A, . 

We consider the words u = 11, v = 01, and the rational numbers X2K+2 and y2K whose binary 
expansions are {0.uv^)2 and (0.11^)2 respectively. According to the functional equation satisfied 
by (Sif ), we have 



S2K+2{^2K+2) = R-Ty/sVo + R-^/3 {R-w/iVo + -R^/3'S'2if (?/2k)) 



-7r/3 ■ 



■Vb- 



■ R2Tv/3S2K{y2K), 



S2K{y2K) = R-w/3 (-R-tt/S^O + RTT/3S2K-2{y2K-2)) 



V_ 



27r/3 ■ 



(2/2 



^2K-2(y2K-2 



2). 



With So(£) = Vb (s is the empty word), we obtain by induction S2K(y2K) = KV_ 



27r/3 ■ 



Vb and 



S2K+2{x2K+2) = {K + l)Vb + y_,/3 + ^2^3 = {K + l)Vb. 

Hence the 0{K\^) we have found in Lemma[7]is satisfying. 

The orbit of the parameterized curve Sk is illustrated in Fig. [2] In this example the recurrences 
which define the sequence (S/f ) and the sequence (F^) are exactly the same. Moreover computing 
the values Sk{x) for a rational x whose binary expansion has length at most K is almost the same 
as computing the values F(x) for a hypothetical solution of the basic dilation equation. This is 
not really true because AoVb is not Vb, and the less significant bits with value must be taken 
into account. (Look at Fig. [2] From the first to the second picture a half-turn has been applied 
to the little pattern of the leftmost picture, because Ag = — I2.) But Ag = I2 and if we use radix 
2^ = 64 in place of 2 this problem disappears. The computation we have made shows that the 



basic dilation equation has no solution: if F was a solution, it would satisfy 'F{y) = V_2,r/3 
for y = 1/3 = (0. 010101... )2. 

We generalize slightly the previous lemma into the following assertion. 



'F{y) 
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Proposition 5. Under Hypothesis |(gev) with < p < A*, the sequence of run 
^,v,^ ^„,^c- fa, A ^o^^^,v,/oW +^ T/(i^-i) o^-A,-Jfic^ IIG„II — n(\K\ f„^ ^„,^„.o, \ \ \ 



nm(; SMms (Si^) associated to V^^ ^satisfies ||Si<-||^ = O(A^) /or every A > A* 
Moreover if X^, is attained, we may replace 0{\^) by O(A^) in case p < A, aKC? &j/ 
0{\fK'') otherwise. 

Proof. The generalized eigenvector \/(''-i) satisfies Q^F^'^^i) = 0{p^K''-^). As 
above we conclude ||Sif ||^ ~ 0{\^) for every A > A*. We may refine the compu- 
tation in case A* is some At : in that case, if p < A* we have HSicHj^^ — O(Af-), 
and if p = A, we obtain ||Sk||^ = OiX^K"). D 

Example 6 [Po^vers of 2 characteristic sequence]. The sequence (n(n)) which takes value 1 
on the powers of 2 and on the others integers admits the generating family («(n), u{2n + 1)) 
and the linear representation 

-o o^ ..=(; ;), * = (; i)^ c=(; 

Matrix Q = Aq + Ai has 1 as a unique eigenvalue and the plane R^ is the characteristic subspace 
Ker(Q — 12)'^ • Hypothesis | ( apt)] is satisfied with 

p=l, R(K) =K, V = { 1 )". 

The joint spectral radius is A* = 1 (consider Aq). The fixed point equation is 

T^)\ 17 )r, \ forO<x<-; <^ „)( ', „ ,„ _,, for - < X < 1. 

F2(x) = F2(2x), - 2 1 F2(a;) = 1 + Fi(2x- 1), 2" 

It admits the only solution (consider F\ and first the interval [1/2, 1[, next the interval [1/4, l/2[, 
and so on) 

(0 1 )*' if <a; < 1, 



^*^^^ \ ( Q f if x = 0. 
(We do not forget the conditions F(0) = 0, F(l) = V.) The sequence (F/f ) satisfies 
(0 0)"' if < i: < 1/2^^' 



Fjf (i) = { ( XjK k/K )*'' if 2'=-V2^' < X' < 2'=/2'f^' for 1 < fc < _ft', 
( l/K 1 )'"■ if z = 1. 

It turns out that (F/^) converges simply towards F. This implies that wc have Sk{x) = K'F{x) + 
0(1) and IISa'IL ~ ■^■ 

Previous Examples [5] and [S] show that under the condition p = A* the sequence 
of running sums {Sk) can or cannot define a limit function F according to the case. 

4.3. General case. The general case is obtained by expanding the vector C over 
a Jordan basis for Q. Let us recall that the height of a generalized eigenvector is 
the dimension of the subspace it spans under the action of Q. Gathering Prop. [4] 
and El we arrive at the following assertion. 

Theorem B. Let L, (^r)o<r<s, C he a linear representation of a rational formal 
power series. The sequence of running sums 

Sk{x) = Y. ^-^ 

\w\ = K 

admits an asymptotic expansion with error term 0{X^) for every A > A*, where A* 
is the joint spectral radius of the family (/lr)o<r<B- The used asymptotic scale is 
the family of sequences p^(^), p > 0, £ (E N>q. The error term is uniform with 
respect to x ^ [0; l]• 
// C is a generalized eigenvector V for Q of height v associated to an eigen- 
value puj of modulus p > A* and if {V^^')o<j<^i, is the associated family of vectors 
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such that V^"-^^ = V, gy(j') = /9wy(j) + V^^-^^ for j > and QV^^'^ = pujV^°^ , 
from the expansion 

is deduced the asymptotic expansion 

\w\=K ^ ^ ^ ^ 

+ • • • + (f ) ipu;f''F<-''~^\x) + (pc^)^r(''-i)(a;) + O(A^) 

where the family (F"')o<j<iy is the unique solution to the system of dilation equa- 
tions 



ri<xi 



ri<a:i 
ri<2:i 

Moreover if the joint spectral radius A* is attained, the error term may be improved 
in the following way. Let V he a Jordan basis for Q. Let m be the maximal height 
of the generalized eigenvectors V from, V associated to an eigenvalue of modulus A* 
such that the vector C has a nonzero component over V , with m — if there does 
not exist such generalized eigenvector. The error term may be taken as 0{\^ K"^). 

The asymptotic expansion is concretely obtained by the algorithm of page 1241 
named LRtoAEl (for linear representation to asymptotic expansion). Clearly this 
gives an asymptotic expansion for the number- valued sum Sk{x) — LSk{x) too. 
It is noteworthy that two important hypotheses are made in the algorithm. First it 
is assumed that the joint spectral radius A* is known and even more it is assumed 
that we know if the joint spectral radius is attained. This is not an obstacle to 
the use of the algorithm. If we have at our disposal an upper bound A-|- > A*, 
we may replace A* by A+, obtaining of course a less precise expansion. Second it 
is assumed that we can solve dilation equations according to Lemma S] This is 
evidently wrong in full generality, but the solution is mathematically well defined 
and may be computed exactly over a dense subset by the cascade algorithm. 

Example 7 [Lipmaa— Wallen's formal power series]. For a nonnegative integer K, the ring 
Z/2^Z may be equipped with two additions. The first one is the quotient addition -|- which 
comes form the addition of the ring of integers Z. The second one uses the identification of 
the set Z/2^Z with the set (Z/2Z)^ through the correspondence which maps a ii'-tuple of bits 
{^k)o<k<K onto the integer eo + 2ei -I- ■ ■ ■ -I- 2^~^eK—l. It is the bitwise addition ©, or exclusive- 
or. The cryptanalysts Lipmaa and Wallen have studied the propagation of differences (in the 
sense of -f ) through © and more precisely the additive difi^erential probabihty adp which maps a 
triple (a, (3, 7) e (Z/2^Z)^ onto the probability that the difference {{x + a) ® {y + fi)) - x ® y 
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Algorithm 1: LRtoAEl 



Input : A linear representation L, {Ar)o<r<B, G , and its joint spectral 

radius A*. 
Output: An asymptotic expansion of the running sum Sk{x) associated to 

the vector C when K goes to infinity at the precision 0~{X^) if A* 

is attained and 0{\^) for a A slightly greater than A* if not, with 

respect to the scale p^ (^g), p > 0, £ e N>o- 

1 Q ■= Z^O<r<B ^r', 

2 if A* is attained then 

3 A '.^ A^ 

4 else 

5 A := any number between A, and the infimum of the modulus of 
eigenvalues of Q greater than A, 

6 end 

7 compute a Jordan basis V for the matrix Q; 

8 expand the column vector C of the linear representation over the Jordan 
basis, as C = J2vev"fvV] 

9 V> := the set of generalized eigenvectors in V such that jv y^ and the 
associated eigenvalue puj has a modulus p > X; 

10 for each vector V in V> do 

11 puj := the eigenvalue associated to V; 

12 ly := the height of V; 

13 an asymptotic expansion 

AvM^) ■■= (/_\pu.)''-^^'F(°Hx)+(f^\pu;f-'^+^F('Hx) + 



■■■+ f^^] {pu;)^-^F^'^-'\x) + (pc^)^F(^-i)(a;) 



is deduced from Proposition [4] 

14 end 

15 Ak{x) := Evev> Jv-Av:k{x); 

16 error termx(a;) := O(A^); 

17 if A, is attained then 
V= :— the set of generalized eigenvectors in V such that 7y / and the 
associated eigenvalue pcu has a modulus p = A, ; 
m := 0; 
for each V in V= do 

u := the height of V; 

m :~ max(TO, v) 
end 
error termi^:(x) :— 0{X^ K™) 



19 
20 
21 
22 
23 
24 

25 end 

26 return 



Sk{x) = A/i" (x) + error termi^: (x) 

K^ + OQ 
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Figure 3. Refinement functions which originate from cryptanal- 
ysis (Ex. [21). 



has value 7 for x, y taken from Z/2^Z with equiprobabihty. By coding eve ry triple (a,/9,7) as a 
word w over the alphabet 'Lj&'L (with w^. = A^a^ + 2/3fc + 7^ for < A; < /^^. iLipmaa et al.l ll2004h 
show that adp turns out to be a rational series. With (_Ej)i<j<8 the canonical basis of Q*, it 
admits the 8-dimensional linear representation 



0<r<8 



(^r)o<r 



C = El 



where Aq is defined as 



Ao 





/ 4 








1 





1 


1 


° ^ 













1 





1 



















1 








1 





1 











1 














4 

















1 


1 






















1 




























1 







V 




















I 



and the other matrices Ar with r ^ are obtained from Aq by permuting row i with row « © r 
and column j with column j r. Using the maximum absolute column sum norm, we find A» < 1 
and because 1 is an eigenvalue for all matrices Ar, < r < 8, we have A, = 1. 
The matrix 



■ Ao+A-i■ 
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is diagonalizable, with eigenvalues 4, 2, 1, and eigenspaces respectively of dimension 1, 3, 4. The 
vector C = El decomposes as C = V4 + V2 + Vi , where 







/ 1 \ 
1 






1 


V4 = 


1 

" 8 


1 
1 
1 
1 
\ 1 / 



V2 





( ' ^ 




1 




1 


1 


-1 


8 


1 




-1 




-1 




^ -3 J 



-=I 



/ 2 \ 
-1 


-1 





-1 






V 1 / 



are respectively eigenvectors for the eigenvalues 4, 2, 1. According to Theorem [B] we obtain the 
following asymptotic expansion (with Fp = LFp for p = 4, 2) 

y adp{«>) = 4'^F4ix) + 2^F2ix) + 0(K), 

(O.TJj)8 <x 

where the functions F4 and F2 are determined by the dilation equation (as usual xi is the first 
octal digits after the octal point of a; £ [0, 1)) 



Fp{x) = - y ArVp + -A^^Fp{8x - xi), with Fp(0) = and Fp(l) 



Vo 



r^x-i 



for p = 4, 2. Functions F4 and F2 are respectively Holder with exponent logg 4 = 2/3 and 
logg 2 = 1/3. Function F^ = LF4 is very near the identity and the graph of a; 1— > F4,{x) — a; is 
drawn on the left of Fig. |3] On the right is the graph of F2 = -LF2. Both functions have an 
amplitude of order about 0.05. The obtained result improves IIDumas et al.L |2007|) which gives 
only an error term 0(2^). 



5. Radix-rational sequences and periodic functions 

It is usual that sequences rational with respect to a radix introduce periodic 
functions in logarithmic scale. These periodic functions appear either by subtle 
elementary arguments, like explicit expression of the sequence with fractional part 
of a logarithm in base B, or by sophisticated argument which express partial sums 
by contour integral of a meromorphic function whose poles are regularly spaced on 
a vertical line. With the approach used in this paper, the introduction of periodic 
functions is natural. Nevertheless these functions are generally only pseudo-periodic 
functions and we will see why they are periodic in the previous works on the subject. 



5.1. Asymptotic expansion. To a rational series S defined by a linear repre- 
sentation L, Aq, . . ., Ab-1, C is associated a sequence u defined in the following 
way: we write the integer n with respect to the radix B, n = {w)b, and the value 
of the sequence for n is the value of the rational series for the word w, that is 
u{n) — (S, w). We want to evaluate the partial sums X]o<n<Af ^(^i) of the sequence 
for N large. With t G [0, 1) and K integer, we may write N = i?^+*, that is K 
is the integer part of log^ N and t is its fractional part. The partial sum up to N 
is the sum of the partial sum up to i?^^ — 1 and a complementary term. The first 
appears as the sum of {S, w) over all words w of length not greater than K minus 
the sum over all the previous words which begin with a 0. The complementary 
term appears as the running sum associated to S we have studied, but without the 
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words of length K + 1 which begin by a 0. This remark gives the following formula 



2<B^+* a<k<K^\w\ = k 



\w' I— fc — 1 



E LA^,C- E LAoA^,C 



\w\=K+l 



\w'\=K 



that is 
(13) 



E uin)^L{U-Ao) E Q'^ 



E 



Ij/i.'inLy . 



i<B^'+' 



0<fe<i<' 



|u)|=A-+l 



Just as we have considered the sums Sk{x), we introduce the vector- valued sums 
(14) 



{w)b<N 



that is 
(15) 



■'N 



iU-Ao) E Q'^- 



E 



AmCy 



0<k<K 



|u;|=_R'+l 



with N = B^^^ , < t < 1. As in the previous section, we decompose the column 
vector C over a Jordan basis for the matrix Q. Hence we consider for the rest 
of the section a family of linearly independent vectors {V'^^^)a<j^i, which satisfies 
gyO") = puV^J) + y(J-i) for j > andQV^) = pujV^, with p > A* and |w| = 1. 
In other words, we assume Hypotheses | (gev)] and (jsr) We want to determine an 
asymptotic expansion for the running sums Sat associated to T^ = v'^^~^\ 

Formula (jlSp expresses the running sums of the sequence as the sum of two 
terms. For the first term in the right member of Formula P^ . we use formulae like 
Eq. (jlOp and the following elementary summation formula 
(16) 

K \ (K'^ 



K 

E 



a 



k~e 



e + 1 

a — 1 ^-^ 



if a = 1, 



K \ a^-'+i 
£ - j) (a - 1)^+1 

1 



+ (-1)^+1 



{a - 1)^+1 



if a^^ 1. 



As a consequence, the first term of Formula (|15[) has an asymptotic expansion in 
the asymptotic scale p^ (,), p> 0, £ ^ N>o- 



28 



PHILIPPE DUMAS 



Lemma 8. Under Hypothesis (gev) with C = V'^'^ ^•' the first term in the right 
member of Formula (|15|) has the following expansion 



K 



(17) (Id-Ao)^Q'=F(^"i) = 



fc=0 

+ E 

€=0 



K \ (puj) 



v-l 



K-u+2 



puj — \ 



-(Id-Ao)F(") 



K 



(puj) 



K-e 



pLJ 



puj — 1 



(Irf-^)F(''-^-i) 






'^-i-i) 



j=i+ 



€=0 ^'^ ' 



in case pu ^ 1 and 

K 



(18) (Id -Ao) 5] Q'^F^^-i) = (Id -Ao)F(°) 



fc=0 



K 



K 



^r; (id-Ao) F('^-^)+F(-^-i +(i,-Ao)F(-i) 



m case puj 



The previous expansion is not an asymptotic expansion but it may easily converted 
into such an expansion in the scale /9^(^), p > 0, i £ N>o, with K = [log^ A^J. 
The error term depends on the relative places of A,, p, and 1. We will return on 
this point in a while. 

The second term in the right member of Formula (fT5)l is nothing but SA'+i(i?*^^)- 
Proposition [4] translates into the following result. 

Lemma 9. Under Hypotheses (gev) and (jsr)| the second term in the right mem- 
ber of Formula (|15p for the running sum associated to V^'^'^^' admits the follow- 
ing asymptotic expansion with respect to the scale p^(^), p > 0, £ G N>Oi with 
K = [log^ iVj , where the functions Y^^' 's are defined in Lemma^ 



(19) 



E 



A^V'^''-'^ 



\w\=K+l 
(u))b<B^+1b'-i 



N- 



K 



(pc^)^-^+2fW(B*-1) 



EfT)(p-) 



(.=0 



K-l 



pujY'^-t-^){B'-^) + f('^-^-2)(B*-1)1 + O(A^). 



for every A > A, and the big oh is uniform with respect to t in [0, 1]. Moreover if 
A, is attained we may replace A by A* in the asymptotic expansion. 

Proof. Proposition [4] and the fundamental formula of binomial coefficients provide 
the expansion of Sa'+i(_B*^^) by a mere substitution. D 

Gathering the results of the previous discussion, we arrive at the following qual- 
itative theorem. 
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Theorem C. Let L, (Ar)o<r<B, C be a linear representation of a radix-rational 
sequence (un). The vector-valued running sum Sat defined by Equation ()14p admits 
an asymptotic expansion with error term 0(7V'°Sb-^) for every A > A*, where A* is 
the joint spectral radius of the family {Ar)o<r<B- If ^* ^s attained the error term 
may be replaced by 0~(Af'°SB '^*), where the soft big oh hides a power of log N. The 
used asymptotic scale is the family of sequences N"'(^'-°^f -I), a G R, £ G N>o. 

Theorem [Cj translates concretely into Algorithm LRtoAE2 of pages [501 [5T] It is 
a simple variation on Algorithm LRtoAEl, but is a little more complicated because 
it takes into account the relative positions of the parameter A, which governs the 
error term, and of the number 1. The right member of Eq. (|19p causes no problem 
because it is readily seen as an asymptotic expansion in the scale N°'['-°^f ^) . The 
right member of Eq. (J17p and (|18p is more annoying. In case A > 1, we take into 
account only the eigenvalues with modulus p > 1 and the case put — 1 oi Lemma [8] 
cannot happen. For the case p > 1, we consider in the right member of Eq. (fT7|) 
only the quantity 



yv — 1/ puj — 1 



u-2 

E 

1=0 



P^ n A ^^Hu-i-i) 
pu 



-(Irf-^O)^^" 






because the quantity 



>^-E(-i)^"^(^:;^(i^-^o)W- 



e=o 

is 0(1) and negligeable in that case. (See lines [T5HH1 and [20H1I1 of Algorithm 
LRtoAE2.) In case A < 1, all terms must be taken into account but we have to 
distinguish between the cases puj ^ 1 and pus = I. (See lines [TMTBl [T7HT51 and [221" 
[23l ) Lastly in case A = 1, the case pcu = 1 does not appear because it is assumed 
p > \, but it may have an effect on the error term if A* is attained. (See lines [2H2l 
and[29H3S) The case pu ^ 1 provide the term X. (See hues [TMH and [20H2ll) The 
term F, which is 0(1), comes automatically into the error term, which is O(A^) 
that is 0(1). 

The asymptotic expansion has essentially the form 

p>\,,e>o ^ ' uj 

where p is the modulus of any eigenvalue of Q, w is a complex number with modu- 
lus 1, and i is an integer related to the maximal size of the Jordan blocks associated 
to pui. As regards $, we change the status of the variable t. It was only an abbre- 
viation for {logg A''}. Now we see it as a real variable and we substitute to it the 
logarithm base B oi N to obtain the expression of the expansion. As a consequence 
the functions $ appear as 1-periodic functions. 
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Algorithm 2: LRtoAE2 (beginning) 



Input : A linear representation L, {Ar)o<r<B, C , and its joint spectral 
radius A». 

Output; An asymptotic expansion of the running sum Sjv associated to the 
vector C when N goes to infinity at the precision 0~{N^°^b ^* ) if \^ 
is attained and 0(Af'°^^ ^) for a A slightly greater than A* if not, 
with respect to the scale 7V"(L'°«| ^J), « e R, ^ e N>o. 

2 if A* is attained then 

3 A !^ A^ 

4 else 

A := any number between A, and the infimum of the modulus of 
eigenvalues of Q greater than A, 

6 end 

7 compute a Jordan basis V for the matrix Q; 

8 expand the column vector C of the linear representation over the Jordan 
basis, as C = Y^vav^vV; 

9 V> := the set of generalized eigenvectors in V such that 7v 7^ and the 
associated eigenvalue puj has a modulus p > A; 

10 for each vector V in V> do 

11 

12 

13 

14 



16 

17 
18 



19 
20 
21 
22 
23 
24 
25 



pLo := the eigenvalue associated to V\ 
V :— the height of V] 
if pw 7^ 1 then 
using Eq. ^ 



X 



v-2 

E 

f=0 



yv — ij puj — 1 



K 



(puj) 



K-e 



puj 



-{U-Ao)V^'' 



-e-i) 



(-1) 



3-e 



pu — 1 






j{ld-A^)V^'''^-^~^ 



Z ■.= X + Y 



else 



using Eq. ([18 



Z 



u-l 



K 



(I, -^)v^(") + E 7 ^^' -^°) ^^"^"'^ + ^^''"'"' + (^"^ -^o)T^('^^^^ 



end 

if A > 1 then 

I ^V;N '■— X 

else 

I Av;Ar := Z 

end 

according to Eq. (I19I1 . a second term is added 



Av: 



V:N 



Av 



V-N 



(pL^)^-''+2r(")(B*-i) 



K 

£=0 



pwF 



{v~t-l) 



(5*-1)+F(''-^-2)(5 



f-2)/r,t-l^ 



26 end 

27 An := I]yev> TvAy^^v; 
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Algorithm 2: LRtoAE2 (end) 



28 error terniAr :— 0{N^°^b ^); 

29 if A, is attained then 
V= := the set of generalized eigenvectors in V such that 7y 7^ and the 
associated eigenvahie puj has a modulus p = A* ; 
m :— 0; 
for each V in V= do 

I' := the height of V; 
m. := max(7Ti, v) 



31 
32 
33 
34 
35 
36 



end 

error terniAr := 0(7V1°8b ^* log^ N) 



37 end 

38 return 



Stv = Atv + error terniAr 



The asymptotic expansion has variable coefficients ( Bourbakilll976l Chapter V) 



The coefficients are taken from the vector space generated by sequences which write 
(u;L'°Sb JVJ(i)(iog5 iV)) '^ith \uj\ = 1 and <& a 1-periodic function. 

Example 8 [Discrepancy of the van der Corput sequence]. iBeiian and Faurg 1 119771 . Il978h 
show that the discrepancy of the van der Corput sequence satisfies the following recursion 

D(l) = 1, D{2n) = D{n), D(2n + 1) = -(D{n) + D{n + 1) + 1) 

and we add D{0) = 0. Let us recall first that the (binary) van der Corput sequence is defined as 
follows: for an integer n, we write its binary expansion (n^-i . ■ ■nino)2; we reverse it and we 
place it after the binary dot. The real number (O.noni . . . rax— 1)2 is the value Xn of the van der 
Corput sequence for the integer ra. Second the discrepancy of the sequence {x„) is 

i^{n, a, 13) 



D{n) = sup 

0<a</3<l 



- (/3 - a) 



where v{n, a, (3) is the number of terms a:^., 1 < fc < ra, which fall in the interval \a, 0). It measures 
the deviation from the uniform distribution for the sequence (a;„) llNiederreiterlll992l). We want 
to ev aluate more precisely the mean value of the sequence (D(ra)) given by IIBeiian and Faurd . 
Il978l . Th. 3) 



1 ^„, , 1 



N 



E^(")., = Tlog2A^ + 0(l). 

The sequence {D{n)) is 2-rational and admits the following linear representation 

/ 1 1/2 \ / 1/2 \ /I 

L = ( 1 1 ) , Aq = I 1/2 , Ai = 1/2 1 , C = 

V 1/2 1 / V 1/2 1 / \ 

with respect to the generating family (D(ra),D(n + 1),1). Numerical computations lead us to 
think that the joint spectral radius of the representation is A* = 1. To prove this result we use 
the Lie algebra g generated by Aq and Ai (with the usual bracket product [A, _B] = AB — BA). 
The derived series of a Lie algebra g is the filtration g = D'^q D D^q D ■ ■ ■ D -D'^g D ■ ■ ■ defined 
inductively by D^g = [g, g], -D*g = [D'°~^g, D'°~^g]. The Lie algebra g is defined to be solvable 
if D*g = {0} for some nonnegative integer k. According to Lie's theorem, if g is a solvable Lie 
algebra over the algebraically closed field of complex numbers th ere exists a change of basis matrix 
which makes triangular all matrices of any representation of g IjGoze and Khakimdianovl . 120001 ) . 
As a consequence if A is a finite set of real matrices which generates a solvable Lie algebra, the 
joint spectral radius of A. is the maximum of the spectral radii p{A) for A in A. Moreover the 
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joint spectral radius is attained with a Euclidean norm w hich writes x i— > (x*'^Px)^ " where P is 
a symmetric definite positive matrix IIBlondel et all |2005|) . 

To apply this result it remains to compute the derived series of g. We use Bq = 2Ao and 
Bi = 2Ai to avoid denominators. Let g be the Lie algebra generated by Bq and Bi . We find at 
the same time a basis of this algebra and its multiplication table, namely 



[, ] 


Bo 


Bi 


B2 


B3 


Bo 





-B2 


-Bs 


~B2 


Bi 


B2 





Ba 


-B2 


B2 


B-i 


-B4 





-Bf, 


B3 


B2 


B2 


-B5 








/ -1 


1 





B2 = [Bi,Bo] = 


1 


1 







\ -1 


1 








' -1 


-1 


B3 = [B2,Bq] = 


1 


1 




-1 


-3 



with 



and B4 = 2_Bo - 2Bi + Bo,, -B5 = -IBo + 2_Bi - 2B3. Hence g admits the basis (Bq, Si, B2, -B3) 
and D^g admits the basis {B2,B3, B4). Next we find that D^g has dimension 1 and is generated 
by ,83 + _B4 = 2Bo — 2Bi + 2,83. As a consequence D^g = {0} and g is solvable. Further the joint 
spectral radius is A* = 1 and it is attained. Because 1 is a simple eigenvalue of Q = Ao + Ai, the 
error term of the asymptotic expansion for the running sum is 0(log A'^). 
Using the basis {Vi,V^, Vj) with 



Vi 



1/2 

-1/2 





the matrix Q takes the Jordan form 






1/2 





2 1 



1/2 

1/2 





: 

The vector C expands as C = Vi + V2 , and because Vi is related to the eigenvalue 1 it may be 
neglected. We find 

SKix) = -2^^ KF°{x) + 2^ F^{x) 
where F" and F^ are solutions of the dilation equations 

F°{x) = -AoF°{2x), for0<a'<l/2, 

F°{x) = -A0V2 + iF°(2a; - 1), for 1/2 < x < 1; 



F^{x) = --F°{x) + -AoF°{2x), 



for < a; < 1/2, 



Fi(x) = --F°{x) + -AoVi + iFl(2a; - 1), for 1/2 < x < 1 

with the boundary conditions FO(0) = 0, F''{1) = V^, Fi(0) = 0, F1(1) = Vj. It is readily seen 
that F''(a;) = XV2, but F-*^ is not explicit. Eventually we arrive at the formula 

Figure |4] show the empirical periodic function and the comparison between the empirical and 
theoretical periodic functions. 

The coefficients of the asymptotic expansion of the running sum S^v inherit the 
properties of the coefficients of the asymptotic expansion of Sk{x). To see this we 
need to be more precise about the way we obtain the asymptotic expansion. Anew 
we use the same formahsm as in Section WA\ 



\ C'~KrTi/j-r>'-rr\Tjr~f tdtt^tu A"^ ; 
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Figure 4. On the left the empirical (red) periodic function of 
Ex. [5] (mean discrepancy of the van der Corput sequence), and 
on the right the comparison between it and the theoretical (blue) 
periodic function computed by the cascade algorithm. 



Corollary 3. Under Hypotheses (gev) and (jsr) the asymptotic expansion of the 
running sum Sjv associated with V^'^~^' has coefficients which are Holder with 
exponent \ogg{p/X) for A > A* and A may be replaced by A* i/ A* is attained. 

Proof. The only point which is to be verified is the continuity at integers, because 
up to this property the Holderian character is evident (the coefficients are expressed 
as combinations of solutions of dilation equations and functions which are smooth 
except perhaps to integers) . As we have seen the expansion of Sjv comes from both 
terms in the right member of Formula P5|) . The first is (I^ ~Aq) J2o<k<K Q'^V^"'^'' 
and we do not want to expand it here, because the asymptotic expansion varies with 
the relative position of p. A,, and 1. Moreover this is not necessary since it depends 
on t only through K. The expansion of the second term is given by Formula (J19p . 
We note that this expression is merely the last column of the matrix F(i?*^^) j5,+^ 
with the notations used in the proof of Lemma 21 Hence we are considering the 
matrix 

A{t) = {U-Ao) Yl Q''V + ¥{B'-^)j;:i+\ 

0<k<K 

where V is the matrix whose columns are the column vectors F^"^ , V^^\ . . ., T/(''~'^) 
as in the proof of Lemma U) The last column of the matrix A(t) is the regular part 
of the expansion of Sat (that is without the error term), and we want to verify that 
1 1— > A(i) is a continuous function at integers. Let Kq be an integer. When t tends 
towards Kq from above, we have K = [ij = Kq and 



lim A(i) = (Id -Aq 

t^Ko 



n<k<Ka 



'V + ¥{l/B)J 



A'o- 



but according to the dilation equation (fT2)l . V{l/B) = AqVJpJ and we find 
fim A(t) = (Id -Ao) J2 Q'^V + AqVJ^J. 



t^Ko 



0<k<Ko 



On the other side when t tends towards Kq from below, we have K = [t\ = Kq — 1 
and 

hm A(f) = (Id -Ao) Y. Q''V + F(l)JJ^" = (Id-Ao) ^ Q'^V + VJI^. 

i^Ko 0<k<Ko-l 0<k<Ko-l 
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0.75 




0.25 



1.25 




0.75 



Figure 5. The change from the function F{x) (left) to the peri- 
odic function $(t) (right) illustrated for the Coquet sequence. 



The difference is 



lim A{t) 

t^Ko 



tKo 



lim A(t) = (Irf ~Ao)Q''''V + {Aa - WVJ^J 



t^Ko 



but Hypothesis (gev) writes QV = VJp^ and we conclude that the difference is 



the null matrix. 



D 



c^^pi 



Theorem [C] deserves some comments. Functions which roughly write <i>(t) = 
*F(i3*~^) have not the self-similar character of the solution F{x) of a dilation 
equation, because we have cut the piece of the function between and l/B. More- 
over it is possible to use a more ordinary scale by expanding ( ^ °^| ^ ) = (^*^ '^^^ ) 
as a polynomial in log^ N with coefficients in i = {log^ TV}. But doing this we hide 
completely the structure of the asymptotic expansion and the dilation equations 
which are behind it. This is certainly the reason why these equations have not been 
perceived before, even if the word "fractal" is frequently employed. 



Example 9 [Coquet sequence] . The Coquet sequence (iCoquetl . 119831) is defined as u(n) = 
(—1)^2(3"). (We recall that S2(n) is the sum of the bits in the binary expansion of n.) It is 
4-rational and it admits (with (u(n),u{4n + 2),u{4n + 3)) as a generating family) the linear 
representation 



1 1 1 

Ao = ( 





10 

Ai = ( 1 -1 





( 1 1 n- 





Aa = ( 1 1 





C- 




1 






1-11 



The matrix 

/ 2 1 1 

(9 = Ao-HAi + A2+A3 = 1 2 -1 
\ 1 -1 2 

has as eigenvalue 3, which is double, and 0, which is simple. The joint spectral radius of the family 
(Ar)o<r<4 is A* = 1, because with the maximum sum column norm we find Aj = 1 and 1 is an 
eigenvalue of each matrix Ar, < r < 4. The column vector C decomposes as C = V3 + Vo, with 

2/3 \ / 1/3 

V3 = ( 1/3 , Vo = -1/3 

1/3 J V -1/3 
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and Vp is an eigenvector for p = 3, 0. 

Theorem [C] applies and we obtain l ICoQueii. [l983l: iFlaiolet et al.lll994h 

E(_l)S2{3n) _ Af'"^*^ 3^~'f*^F(4'(*''~-'-) + 0(1), 
JV^ + oo 
n<N 

with t = log^N. The function F is the sum F = Fi + F2 + F3, where F = (^1,^2,^3) is the 
unique solution of the dilation equation, written with Convention |(hec)[ 

Fi(x) = -Fi(4x) + -F2(4x) + -F:i{4x) + -Fi(4x - 1), 

F2{x) = -F2{4.x - 1) - -F3(4x - 1) + -Fi(4x - 2) + -F2(4x - 2), 

^3(^) = ^Fsiix - 2) + ^Fi{4x - 3) - ^F2{ix - 3) + ^^3(4^ - 3), 

with the conditions Fi(0) = ^2(0) = ^3(0) = 0, Fi{l) = 2/3, ^2(1) = ^3(1) = 1/3. The func- 
tion F{x) and the periodic function #(t) = 3^~i'j -F(4i*J~^) are illustrated in Fig. [5] re spectively 
on th e left and on the right. The positive character of $ proves a subtle phenomenon IINewmarJ . 
Il969h : the number of ones in the binary expansion of the integers which are multiple of 3 is more 
often even than odd. Both functions are Holder with exponent log4 3 ~ 0.795. In the translation 
from F to <E>, the first quarter of F is lost and <I> has not the auto-similar character of F. 

As we have seen if the previous sections, it may be interesting to change the 
radix. Going from radix B to radix B^ , we change the matrix Q into its power Q^ . 
Frequently the eigenvalues of Q have arguments which are commensurable with tt 
and a suitable choice of T gives a matrix Q whose all eigenvalues are real nonnega- 
tive. This trick is often used implicitly. For example the Coquet sequence i s viewed 
as a 4-rational sequence (as we have made in the prev i ous e xample) in ( Coquetl . 



19831 [Pumont and Thomas. Il989t lAllouche and Shallitl . l2003l ) (but this is not the 



case in (JFlaiolet et al.l . ll994l )'). However it may be of interest to use two radices at 
a time. Even if both dilation equations define as well the function under consider- 
ation, it may be practically simpler to study the solution for one radix in terms of 
solutions for the other radix. The next example illustrates this point. 

Example 10 [Rudin-Shapiro sequence continued]. As an application of Theorem ICl we 
obtain the asymptotic expansion for the Rudin-Shapiro sequence defined in Ex.lTl llBrillhart et al.L 
Il983h 

y n„ = v^-J'(log4Af)-f 0(1) 

where $ is the 1-periodic function defined by #(i) = 2^~'*'_F(4i'j ^) and F is defined through 
a dilation equation for radix 4. Functions F and <& are illustrated in Fig. \E\ Let us denote ttj; 
the part of the graph of F which corresponds to the interval [fc/8, (k + l)/8] for < fc < 8. It 
is evident that the parts of odd index on one side and the part of even index on the other side 
reproduce the same pattern (with a piece upside down) and we want to prove these facts. 

Besides of radix 4, we will use radix 2. The matrix Q for radix 2 (see Ex. [TJ has two eigenvalues 
it V2 and associated eigenvectors V± = (2it \/2, it v2)/4. The vector C of the representation writes 
C = V+ + V— and we have to consider two vector- valued functions F±. Their components are -F±,i 
and F±^2- The function under study is _F = F+,i + -F-1-,2 + ^-,1 + ^-,2. 

To abbreviate the computations and with the hope to make them clearer, we introduce the 
following formalism. The group of affine transforms of the real line acts on the vector space of 
functions from the real line into itself. More precisely the group acts on the right by substitution: 
for an affine transform a and a function /, the image f.a is simply the composed map f o a. As 
affine transforms it is natural in the actual context to consider do : x 1-^ 2x and di : a:: 1— > 2a; — 1. 
We will use ti/2 : x 1— > x -t 1/2 too. With these notations the dilation equation for F± writes 
(We use Convention | (hec)] to write the dilation equations, that is we consider that functions are 
adequately extended to the whole real line.) 

(20) F±,i = it-pF±,i.doit-pF±,2.rfo, F±,2 = it-pF+,i.diT-pi^±,2.'ii, 
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Figure 6. The solution of the dilation equation (left-hand side) 
and the periodic function (right-hand side) associated to the 
Rudin-Shapiro sequence show some symmetries. 



with the conditions F±(0) = and F±(l) = V± We will collect some simple facts in order to 
achieve our goal. It results from the dilation equations that F±,2 is on the left of 1/2, while 
F±^i is constant on the right of 1/2. The action of tx/2 o^i both members of the second equation 
above with the equality di o tj ^ = d,Q provides 

A subtraction and an addition give respectively 

(21) F+,1 - F+,2.ti/2 = V2F+^2-do, F+,1 + F+^2-h/2 = V2F+,i.do. 

From these preliminaries follow the following formulae 



(22) 



F(a;+i) - F(x) 
F{x+^)-Fix) 



for < X < J , 



for < X < 



F {x + ^) + F{x) = 2 
F{x+\)+Fix) = l 



for i < X < i, 

for — < a:: < — . 

8 ~ ~ 4 



We prove only the third one; the other proof are of the same type. For < a; < 1/8, the number 
X + 1/4 lies in [1/4, 3/8] and both numbers are in [0, 1/2]. The dilation equations I I20I I give 



F±,i(x) = ±^F±,ii2x) ± ^F±,2(2x), 



F±,2{x) = 0, 



and 



F±,i X 



■^7^2^^- 



('^ + D^71^^'^''"- 



F±.2 X 



Since F is F+,i -I- i^+,2 -I- -F'-.l -I- -^-,2, we obtain 



F { X 



F{X) : 



2x- 



V2 



F+,1 2x 



-F+'2 \2x + 



F_,2 2x + 



F+^{2x) - F+,2(2x) + F_,i(2a;) + F-^i'^x) 



However 2x is in [0, 1/4] and 2a; -1-1/2 is in [1/2, 3/4], The functions F±,l are constant on the right 
of 1/2, while the functions -F±,2 are constant on the left of 1/2. Making explicit theses constant 
values, we obtain 



F X 



■ ^{=^0 



1 

72 



-V2- F+^x(2x) + F_,i(2x) -I- F+ 



2X' 



2x- 



ASYMPTOTIC BEHAVIOUR OF RADIX-RATIONAL SEQUENCES 37 

Besides Formula l|21|l gives 

F±,i{2x) = F±,2 h^+l 

because x is in [0, 1/8]. The expected formula is proved. 

The previous computations ask some questions. We have established formulae which write 

f.a 

where / runs through the four functions i^+,i, F+,2, F—,i, F-,2\ o, runs through the set of affine 
transforms of the real line; the family (Aj „) has a finite support; c is some constant. Is it possible 
to design a method in order to find such relationships? For example, how to obtain mechanically 
the formula 

^F±a(.)--lF,,.(l-.) = 17 
Does exists a regular process which gives formulae of the type 

Y^ XaF.a = C? 

a 

Evidently we want a minimal generating set for such formula;. 

From the four formulae II22I I above we deduce that all pieces of the graph of F are obtained 
from TTQ and tti by some translations or glide reflections. The symmetries of F are translated 
to <&, but they lose their graphical evidence. The pieces ttq and tti disappear and the pieces nf;, 
2 < k < 8 become the pieces 7r^ associated to the interval [log4(A:/2),log4((A; + l)/2)]. The links 
between the pieces become more intricate. For example the pieces tTj and TTg on one side and TTg 
and TT^ on the other side are linked by the formula 2''<l>(s) + 2'<I>(t) = 4 under the condition that 
both numbers s e [0, 1/2] and t e [log4 3, 1] are related by 4' - 4'' = 2. 

Example 11 [Rescaling]. In the study of radix-rational sequences, it is an attractive idea at 
first sight to multiply all the matrices Aq, . . ., Ab—i of a linear representation by a scalar in order 
to control the order of growth of the sequence. But after a while the idea seems to be silly, because 
it introduces a weighting according to the length of the radix B expansions of the integers. The 
previous theorem shows clearly the effect of such a rescaling. It appears that the modification 
does not change at all the functions F{x) defined by a dilation equation, and changes very simply 
the periodic functions $(i). 

Let us consider the sequence u{n) = n. It admits the representation 

L=(0 1), Ao=(f ;), ^,. = (f ;)for.>0. C=(; 

for radix B. There is a dominant eigenvalue p = B^ , with associated functions F{x) = x^ and 
^(t) = 1/2. We obtain the obvious asymptotic expansion 



y n = -N'^ + 0(N). 



N^ + oo 2 
n<iV 

Now consider the sequence u{n) = n/B^B (") ^ where A^ (n) is the length of the radix B expansion 
of n. With S = 2 it begins with 0, 0.1, 0.10, 0.11, 0.100, 0.101, 0.110, 0.111, 0.1000, ... It admits 
the representation 

L=(0 1), Ao=(; J^), A.= {,.jj, ,;^)forr>0, C = (^ ^ 

and again we have F{x) = x"^ , but p = B. We obtain the less obvious asymptotic expansion 

y ", , = -Af<I>(logs Af) + 0(1) 

n<N 

with 

<!>(*) = = — — - cosh (t ) In B , for < t < 1. 

The convergence towards this function is illustrated in Fig. [7] where a (piece of) catenary appears 
in an evident way. 
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Figure 7. The empirical periodic function for sequence u{n) = 
j^^2'^2(n) (where X2{n) = 1 + [fog2 nj is the length of the binary 
expansion of n) as a function of logj n. A catenary appears. 



5.2. Substitutions and automata. In the very beginning of the introduction 
we have written that nobody has asserted a general theorem about the asymp- 
totic expansion of radix-ra tional sequences. This is not quite the truth because 
AUouche and Shallid ([20031, Sec. 3.5) provide a theorem. But their statement suf- 
fers from a severe restriction because they assume for example the hypothesis 
Aw = 0{K) for \'w\ = K, that is A* < 1, and some other technical hypotheses. 
It is noteworthy that the res ult they give ha s no error term and for this it appears 



as a direct generalization of (|Delangd . 119751 ) 



Other work s which deserve attention ar e related to the name of Dumont ( Dumont and ThomasL 
19891 : iDumontl . Il990t iDumont et al.l . Il999r ) . There is a link between these works and 
the present one, but it needs some explanations to be emphasized. We consider a 
finite alphabet X. A substitution ct on ^ is a map which associates to each letter 
a non empty word. It extends as a morphism of X* by concatenation. If a is a 
letter which turns out to be a strict prefix of the word a {a), the sequence (cr'^(a)), 
obtained by iteration of a fr om a, conver ges to an infinite word co — {i-Un)n>Oi which 
is named a fixed point of a (|Fogeil2002f) . 

We may associate to the substitution a and to the letter a an automaton Aa- 
The state of the automaton are the letters of the alphabet. There is a unique initial 
state which is the letter a. All states are final states. The transitions are labeled 
by nonnegative integers and obtained in the following way. For a letter x, the 
non empty word a{x) has length £ > 1 and writes a{x) = ao{x)ai{x) . . . ai-i{x). 
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There is a transition from state x to state ar{x) labeled by r. Conversely the 
automaton determines the substitution and the fixed point of the substitution which 
begins by a. If B is the maximal length of the words <7{x) for x a letter, all labels 
are bounded above by i? — 1. To the automaton is associated the language C 
it recognizes, that is the set of words over the alphabet {0,1, . . . ,B — 1} which 
translate a sequence of transitions from the initial state to a final state. This 
language is regular because it is recognized by the automaton. The simplest case is 
the case where the substitution is of constant length, which means that all words 
<j{x), X € X, have the same length B >2. In such a case the language C is merely 
the free monoid {0,1, . . . ,B — 1}*. 

The automaton may be enriched by an output function / which associates to 
each state, that is letter, a value. The output function ext ends to words and the 
i mage of the fixed point lu is an infinite sequence of values. iDumont and Thomad 
( 19891 ) study the asymptotic behaviour of the running sum 



S^N)^ ^ /K) 



0<n<Af 



where / is real-valued, under s ome hypotheses which we wi ll describe later. 

To merge the framework of ( Dumont and Thomad . Il989l ) and ours it suffices to 



consider a linear version of the automaton. For a given commutative field K, we 
consider the space K'*, whose canonical base is {ex)x£X- To each integer r between 
and i? — 1 we associate a, d x d matrix Ar, where d is the size of the alphabet, 
defined by Arex = Ry if there is a transition labeled r from state x to state y in 
the automaton, and A^Cx = otherwise. The matrix Q = A^ + ■ ■ ■ -\- Ab-i is the 
incidence matrix of the substitution; its entries Qx,y is the number of occurrences of 
the letter y in the word a{x). We insert in the description the column vector C = Ca 
and the row vector L = {f{x))x<£A and we have a linear representation. For a word 
w — wq ■ ■ ■ wk-1 taken from the language C of the automaton, it is equivalent to 
follow from the initial state a the transitions labeled wo, ■ ■ ■, wk-i and to compute 
the value f{z) of the state z where ends the path, or to compute the matrix product 
i^ujif-i • • • A^^C . We are not far to the definition of a radix-rational sequence, but 
the idea of a numeration system is yet lacking. 

The lengths of the iterates cr" (a) are the elements of an increasing sequence B — 
{Bn)n>o which begins with Bq — 1. This sequence defines a numeration system: 
a nonnegative integer n admits an expansion (qk-i ■ ■ ■<Io)b in the system B if it 
writes n = qoBq + aiBi + • ■ • + uk-iBk-i- There is always at least one expansion 
named the normal 6-representation, which is obtained by a greedy algorithm. The 
language S associated to the system of numeration is the set of all normal B- 
representations of nonnegative integers. (The representation of is the empty 
word.) Generally speaking, the language £ is not regular and is distinct form the 
language S, but in the simple case where the substitution is of constant length B, 
both languages C and S are equal. As a consequence the sequence (/(wn)) is 
S-rational. It is even B-automatic: for an nonnegative integer n we write its 
radix B expansion and we use this word to follow a path which begins with a in 
the automaton; the end of the path is w„ and the value associated to n is /(wn). 
It is proved that a B-automatic sequence is a B-rational sequence which takes a 
finite number of values. 

Nevertheless the previous system of numeration is not suitable in case o f a sub- 
stitution which is not of constant length. The prefix-suffix automaton A'^^ (JMossa . 
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19961 : ICanterini and Siegej l200lh is the right automaton. Its states are the letters 



of X. There is a transition labeled (p, x, s) G X* x X x X* from x to y if cr(y) — pxs 
(hence the name of the automaton). All states are initial and a is the only final 
states. Clearly A'^ is a sophisticated version of A (the transitions are reversed but 
this is of no importance), whi c h is e quivalent in case of a substitution of constant 
length. Dumont and Thomas (l989j use a version of this automaton which takes 
into account only the prefixes. The language C recognized by A!^ may be defined by 
a set of prohibited patterns: a pair of labels ((p, x, s), (jp' , y, s')) such that a{y) / psx 
does not occur in a word of C . For every subword uj^ — {oL!n)o<n<N there exists 
a unique path {{pk,Xk,Sk))o<k<K of length K, determined by Bk < N < Bk+i, 
which ends at a. In other terms there exists a unique word {{pk,Xk-, Sk))Q<k<K 
from C of length K such that Xk-i = a- It defines a cutting of uj-'^ = {uJn)o<n<N 
which writes uj^ = a^~^{pK-i)<y^~^{pK-2) ■ ■ ■O'°(po)- Moreover the prefixes pk 
lie in the finite set of strict prefixe s of the words (j{x) with x a letter. Such a result 
demands hypotheses ( Mossall992i ): the substitution is primitive, which means that 



there exists a k such that for every pair of letters x and y the letter x occurs in 
(T (y); the fixed point uj is not a periodic sequence. Taking into account the lengths 

of the words, wc obtain TV = \a^-'^{pK-i)\ + \(7^~^{pK-2)\^ \-\a°{po)\- This is 

the expansion of the integer N according to the numeration system associated with 
the substitution a and a. (Note that there does not exist a general definition of 
what is a numeratio n system.) The proof of the existence of the expansion is anew 
a greedy algorithm ( Dumont and Thomad . ll989L Lemma 1.3): px-i is the longest 



prefix of oj such that a^~^{pK-i) is a prefix of co^ . 

The asymptotic study of the sequence {S^{N))n>o is based on the following 
remarks. The function / is defined on the letters of the alphabet X and extends 
additively to the word of X* . Particularly we have S^ {N) — /{uj^). Let L the row 
vector whose entries are the values of / on the letters of the alphabet; let Cw the 
column whose entries are the number of occurrences of each letter in the word w. A 
consequence of these definition is the formula f{w) — LCw Let Q be the incidence 
matrix of the substitution a. We see immediately the recursion C„i,jj\ — QCw As 
a consequence we have two analogous formulae 

K-l K-l 

fe=o k=a 

where U is the row vector w hose all entries are equal to 1. Both formulas render 
plausible the following result ( Dumont and Thomad . Il989l . Th. 2.6), 



S^{N) = LAN + \ogg NN^G{N) + o{logg NNf^). 

It is proved under the following hypotheses and through the following assertions. 
Because Q is a primitive nonnegative matrix, its spectral radius is a dominant 
eigenvalue. The column vector A is its positive eigenvector, normalized by C/A = 1. 
The existence of a sub-dominant eigenvalue 0' is assumed. It dominates all others 
eigenvalues apart 9, that is > 9' > \9"\ for all eigenvalues 9" except 9 and 9'. 
Moreover it is assumed 9' > 1. The integer a is defined such that a + 1 is the 
multiplicity of 9' as root of the minimal polynomial of Q and the real number (3 is 
defined by /3 — \ogg 9' . The function G is defined as the limit 

_.-/3 K^ Sf{\_9'^x\)-LK\9^x\ 



Vx > 0, G{x) = X-'-' hm 



fc^+oo k°'9'^ 
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It satisfies G(6x ) — G(x) for x positive and is Holder of exponent f3. 
To summarize lDumont and Thomad ( 19891 ) study the mean asymptotic behaviour 



of real- valued sequences defined by substitutions and for substitutions with constant 
length this reduces to automatic sequences. We study complex-valued sequences 
associated to complex rational formal power series and for series which takes only 
a finite number of values this reduces to automatic sequences. 

5.3. Linear representation insensitive to the leftmost zeroes. In the use of 

Formula (J13p . we may introduce a simplifying hypothesis which turns out to be non 
restrictive. To a rational formal series we associate a radix-rational series through 
a numeration system. But the knowledge of a radix-rational sequence does not 
determine a unique rational formal series, because the expansions of the integers 
does not use the words which begin with a zero. A first idea which comes in mind 
to complete the definition of the formal series is to decide that words which begin 
with zero give a null value. But there is another way to extend a radix-rational 
sequence into a formal series and it is more natural after all. Let us say that a 
linear representation of a radix-rational series is insensitive to the leftmost zeroes if 
it satisfies LAq = L. For such a representation the formal power series is completely 
determined by the radix-rational series, because the adding of some zeroes on the 
left of the expansion of an integer does not change the value associated to this word. 

Lemma 10. Every radix-rational sequence has a (reduced) linear insensitive to the 
leftmost zeroes representation, that is a representation such that LAq — L. 

This point is described in f Pumad . Il993l sec. 4.2) (where an insensitive to the 
leftmost zeroes representation is termed standard) , but we give a proof because it 
explains how most of the linear representations in this paper are obtained. 

Proof. For the sake of simplicity, let us assume that the radix is B = 2. The 
hypothesis is that all subsequences (M2'=n+r)i fc > 0, < r < 2'^, of the radix- 
rational sequence (m„) generate a finite dimensional vector space U. Let d be the 
dimension of that vector space. (For the null sequence, we have d — 0. We exclude 
this case.) We consider the subsequences (u„), (u2„), {u2n+i), (w4„), (u^n+i), 
{uin+2), ■ • . in that order and the dimension dj of the vector space generated by 
the j first subsequences (with d_i = 0). We select the subsequences w^, v^, . . ., 
v'^ such that the dimension increases by 1 when they encountered in the previous 
list. The family {v^)i<i<d is a basis of U. The sequence (w„) expresses as a 
linear combination of w^, w^, . . ., w'' (the sequence (u„) is nothing but v^), and 
this gives the column vector C of the linear representation we are building. Next 
we consider the action of and 1 over the sequences v^ defined by (0.w^)„ = ?;2„ 
and {l.v^)n = ''^2n+i ^^'^ ^^ express the images in the basis (f^)i<^<d. This gives 
the square matrices Aq and Ai of the representation. At last the row vector L 
is the vector of initial values Vq. In that way, we have a linear representation 
of the sequence (««). This may be readily verified by considering the sequence 
of row vectors A„ = (u^, w^, ... u^) which satisfies m„ = A„C A2„ = A„ylo, 
A2n-i-i = A„Ai, L = Aq. (Note that lAUouche and ShallitI (|2003h use a column 



vector for A„, hence the transposed matrices for Aq, Ai.) Moreover the equation 
(O.'(;^)o = Vq, that is AqAq = Aq, renders evident the formula LAq = L. The proof 
would be more enlightening by considering a binary tree and a prefix part of the 
monoid of words, but we will not insist on this point. D 
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Figure 8. The empirical rosette of Ex. [THwith -d = 27r/5. For 
each frond of the rosette, the index N goes from 2^ to 2^+^ — 1. 
Here K goes from 11 to 15 and the color from red to blue. The 
vertices of the pentagon are the points T{K) with K integer. 



We do not have used this hypothesis of insensitivity in the proofs because it 
is unsuitable for the vector-valued version we have elaborated, but most of the 
examples we have given have this property and this saves us to compute the first 
term of Formula (1131) . 



5.4. Periodicity versus pseudo-periodicity. The occurrence of periodic func- 
tions in logarithmic scale is common in the asymptot i c study of radix-rationa l 
sequ e nces (Coquet, 1983; Delangc, 1975; Dumont et al.l . Il999t iFlaiolet and Golin , 
11994 iFlaiolet et all . 11994 lOsbaldestin and Shiul . ll989f) . Nevertheless the rising of 
periodic functions, if true in practical examples, is wrong in full generality. The 
reason is the following: in common examples the eigenvalues of the matrix Q have 
arguments which are commensurable with n. As we have already explained, a 
change of radix permits to consider that the useful eigenvalues are all positive 
numbers. As a consequence functions $ are periodic, because they depend only on 
the fractional part of the variable t. But in the general case they write something 



like $(i) 



;Wpl- 



-{t}p(^B{t} 1) and if u does not write u) = e*'' with ^/tt a ra- 



tional number, the function $ is not periodic but only pseudo-periodic. Below we 
exemplify this phenomenon. 



Ao- 



Example 12 [Rosettes]. Let us consider the linear representation for radix B = 2 and dimen- 
sion 2, 

cosi9 ^ A - f ^ -sini? 

cos)9 J ' ^ ~ \ sini? 

where i? is a real number restricted by the condition i9 ^ mod it/2 to avoid degeneracy. We do 
not define the row vector L, because we are interested in the vector- valued sequence 

Sjv = 2^ Au]C 

0<n<N 

(w is the binary expansion of the integer n). With regard to the column vector C, we note that 
both matrices Aq and A\ write Ao = cosi9 I2, A\ = sini?iJ^/2 (notation of Ex. [5] Rtp is the 
rotation matrix with angle </?). As a consequence they commute with all rotation matrices and 
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if C is changed into R^C, then Sjv is changed into R^Yli^i. Evidently we have a similar formula 
with dilations and so we may take as vector C the first vector Ei of the canonical basis. 

The joint spectral radius is A* = max(| cos»9|, | sini?|), because ||Ao||i = | cosi?|, ||j4i||-^ = | sini9| 
and j4g and A^ admit respectively the eigenvalues cos^ i? and — sin'^ i?. The matrix Q = Aq + Ai 
is the rotation matrix _R^ and its eigenvalues are pu)± with p = 1 and liJ± = e^ . We take as 
eigenvectors respectively for liJ_|_ and tj_ the vectors V^ = ( 1 ~* ) > ^— ~ ( -'^ * ) ■ "^^^ 
vector C expands as C = (V+ + V—)/2. We have A* < p and Theorem [C] applies. The asymptotic 
expansion for the sum S_|_ jv associated with V+ writes 

K 
S+,^r = (l2 -Ao) Y, Q'=y+ + e»('f^'+i)''F+(2{*}-i) + 0(Af ), 

*,— 
where Y± is defined by the dilation equation 



:cosi9eT"'F±,i(2x), 
:cosi?eT'''F±^2(2x), 



1), 



for 



1 



< !• < 1, 



F±,2(3:) = cos,?eT'''F±,2(2x), lor u ;, a; <, ^, 

: COS i?eT*' 
: =picosi9eT' 

with the boundary conditions F±(0) = 0, F±(l) = V±. Both vector- valued functions F± are 
conjugate. Emphasizing their real and imaginary parts U and itV and carrying on with the 
computation we obtain for the sum Sjv associated with C = E\ the asymptotic expansion 



-sin)?e=F»''F±,2(2x- 
F''' -1- sini?eT»'>F±^i(2a; - 1), 



Sat = r(i) + 0(Af ), 

iV— i-l-oo 



with r : t I 



Obviously an increasing of t (and consequently of K) by 1 rotates the vector A(4) by i?. The 
change for B is less obvious. From the real point of view, the dilation equation for F4- rewrites 



A(t)+B(t), 




. -a / sin-d/2 \ ,,, .-Of cosUK + l/2){) ~ 
"" 2 V cos J/2 ) ' ^(*) ="■^2 1 sin({V + 1/2)^ - 


V2) 


B(t) = cos({_ft' + l)i9) U(2<*>-^) - sin((X + 1)0) V(2(*> 


-^). 



#(x') = To#(2x) for < x < 



1 



*(X) : 



where ^ is the 4-dimensional vector-valued function # 
have the following expressions 

3s2i9 



ToW-f Tl#(2a•- 
( U V )' 



1) 



r 1 ^ 

tor — < X < 1, 

2 - 



To: 



V 





- cos i9 sin i9 




cos^ -d 



- cos i9 sin i9 



cos 1? sin 1? 



cos^ i9 





\ 

cos i9 sin 1? 



cos^ i9 



The matrices To and T\ 



-■ cos "d R_ 



Ti 





cos i9 sin i9 



— sin^ 1? 



' cos -d sin -d 



sin^ i9 







sin^ 1? 



cos 1? sin -d 



— sin it 



- cos i9 sin i9 





:sini?i?_ 



)R. 



The boundary conditions are #(0) : 
the matrix 



and *(1) 
( 



W with W -- 
1 



(1 -1 ) 



7r/2- 



Let P be 



-1 \ 
10 
10 

V -1 y 

We verify P-^ToP = To, P'^TiP = Ti, PW = W. As a result P* is a solution of the dilation 
equation. But the dilation equation has a unique solution and we conclude that ^ satisfies 
P# = #. This means —V = R^^2^- This shows the formula B(t -|- 1) = _R^B(t). Eventually we 
see that T{t + 1) is the image of T{t) by the rotation of angle i9 about Q. 

Let us assume that 1? is commensurable with tt, and let us write "S/it = p/q where p/q is in 
lowest terms. We see immediately that the arc T is 2g-periodic, and even 2^~''5-periodic if k is 
the dyadic valuation of p. This phenomenon is illustrated by Fig. |8] Besides it is easy to verify 
r(i-|-2-'*5)-|-r(t) = 2C. 
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To the contrary when the ratio iJ/tt is irrational, the rotation i?^ leaves invariant the image 
of the arc F but its order is infinite. The non-periodic character of the arc T is evident when we 
consider its value at integers, 

T./.^N ^ '^ f cos(K ~ 1/2)1? 

^ ' 2 V sm(K — 1/2)1? 

Hence we see that the functions <& are not necessarily periodic, even if they are frequently periodic 
in concrete examples. 

Acknowledgment. I am most grateful towards Frederic Chyzak (Algorithms 
Project, INRIA) for his ever patient ear during my effort to understand radix- 
rational sequences. 
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